(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     98614,       2088]
NotebookOptionsPosition[     97804,       2057]
NotebookOutlinePosition[     98207,       2074]
CellTagsIndexPosition[     98164,       2071]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"<<", "ACPackages`"}]], "Input",
 CellChangeTimes->{{3.50692072647425*^9, 3.506920728864875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"$WD", "=", "\"\<C:\>\""}], ";"}]], "Input",
 CellChangeTimes->{{3.505358892078125*^9, 3.505358909046875*^9}}],

Cell[BoxData[
 RowBox[{"$WD", "=", 
  RowBox[{"SetWorkDir", "[", "\"\<2014.11.19\\\\2\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.506147355829625*^9, 3.5061473611265*^9}, {
   3.50726751628675*^9, 3.507267517146125*^9}, {3.5376857609796247`*^9, 
   3.5376857676202497`*^9}, 3.538129311765625*^9, {3.599812956215828*^9, 
   3.599812962597664*^9}, 3.5998159499941483`*^9, {3.604125951876953*^9, 
   3.604125954845703*^9}, 3.604282518145508*^9, {3.60430203384375*^9, 
   3.6043020376376953`*^9}, {3.605260018712188*^9, 3.605260021061797*^9}, {
   3.6176233790800624`*^9, 3.6176233858281093`*^9}, {3.623482089254883*^9, 
   3.6234820893876953`*^9}, {3.623484763916992*^9, 3.623484774814453*^9}, {
   3.623487015705078*^9, 3.623487015830078*^9}, {3.624259469510742*^9, 
   3.6242594699746094`*^9}, {3.6246875306953125`*^9, 3.624687535524414*^9}, {
   3.6246976584677734`*^9, 3.6246976621464844`*^9}, {3.6248644395751953`*^9, 
   3.624864439881836*^9}, {3.6248664769179688`*^9, 3.6248664787246094`*^9}, 
   3.6259067990078125`*^9, {3.6260742754179688`*^9, 3.6260742783466797`*^9}, {
   3.626079990010742*^9, 3.6260799902890625`*^9}, {3.6266058525595703`*^9, 
   3.626605857205078*^9}, {3.6266123118066406`*^9, 3.6266123181816406`*^9}, {
   3.6266326809853516`*^9, 3.626632684736328*^9}}],

Cell[BoxData[
 RowBox[{"FileNames", "[", "\"\<energy*\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.6042131756972656`*^9, 3.6042131833242188`*^9}, {
  3.6260742805820312`*^9, 3.6260742813466797`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"SetDirectory", "[", "$WD", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"infos", "=", 
   RowBox[{"SortBy", "[", 
    RowBox[{"ScanSnaps", ",", 
     RowBox[{
      RowBox[{"#", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}], 
      "&"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Prepend", "[", 
   RowBox[{"infos", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<Time\>\"", ",", "\"\<Iteration\>\"", ",", "\"\<Reynolds\>\""}], 
     "}"}]}], "]"}], "//", "TableForm"}]}], "Input",
 CellChangeTimes->{{3.536049279181506*^9, 3.536049286181506*^9}, {
   3.536049317431506*^9, 3.536049328775256*^9}, {3.5376857784483747`*^9, 
   3.5376857849014997`*^9}, {3.5376862506202497`*^9, 
   3.5376862695264997`*^9}, {3.538129304109375*^9, 3.5381293043125*^9}, {
   3.53829224075*^9, 3.5382922416875*^9}, 3.538389701015625*^9, {
   3.538738063171875*^9, 3.538738063390625*^9}, {3.5998129667392654`*^9, 
   3.599812970460945*^9}, {3.5998159523330154`*^9, 3.5998159523906326`*^9}, {
   3.5998949257216873`*^9, 3.5998949269199295`*^9}, {3.6042825645859375`*^9, 
   3.6042825681123047`*^9}, {3.604302048413086*^9, 3.6043020692509766`*^9}}],

Cell[CellGroupData[{

Cell["\:041d\:0435\:0441\:043a\:043e\:043b\:044c\:043a\:043e Re \:043f\:0440\
\:0438 \:043e\:0434\:043d\:043e\:043c \[Kappa]", "Title",
 CellChangeTimes->{{3.5376859382139997`*^9, 3.5376859455108747`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"CellPrint", "[", 
   RowBox[{"Cell", "[", 
    RowBox[{
    "\"\<\:041c\:043e\:0436\:043d\:043e \:0441\:0442\:0438\:0440\:0430\:0442\
\:044c :) \>\"", ",", "\"\<Title\>\""}], "]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\:0420\:0430\:0441\:0441\:043c\:0430\:0442\:0440\:0438\:0432\:0430\
\:0435\:0442\:0441\:044f \[Kappa]=\>\"", ",", 
    RowBox[{"kappa", "=", "0.975"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"len1", "=", 
   RowBox[{"Import", "[", 
    RowBox[{"\"\<energy.dat\>\"", ",", "\"\<Table\>\""}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\:0420\:0430\:043c\:043a\:0438 \:0432\:0440\:0435\:043c\:0435\:043d\
\:0438: \>\"", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"len1", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", 
      RowBox[{"len1", "[", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "1"}], "]"}], "]"}]}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ReadTextSnapFile2D", "[", 
   RowBox[{"\"\<pumpflowRe80.snp\>\"", ",", "5", ",", 
    RowBox[{"PrintInfo", "\[Rule]", "None"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"coor", ",", "bounds", ",", "node"}], "}"}], "=", 
   RowBox[{"ReadAuxFiles2D", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"\"\<coord\>\"", ",", "\"\<node\>\""}], "}"}], ",", 
     RowBox[{"ShiftVector", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}]}], ",", 
     RowBox[{"ProcDistrib", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"8", ",", "16"}], "}"}]}], ",", 
     RowBox[{"NumberOfRefs", "\[Rule]", "2"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"M", "=", 
    RowBox[{"1", "+", 
     RowBox[{"kappa", " ", 
      RowBox[{
      "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}]}]}], 
   ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"rz", "=", 
     RowBox[{"Outer", "[", 
      RowBox[{"List", ",", 
       RowBox[{"coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}],
        ",", 
       RowBox[{
       "coor", "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}]}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"dx", "=", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", 
         RowBox[{"Subtract", "@@", "#"}]}], "&"}], "/@", "bounds"}], 
      RowBox[{"Dimensions", "[", "vx", "]"}]]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Theta]", "=", 
     RowBox[{"Map", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"ArcTan", "[", 
         RowBox[{
          RowBox[{"#", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}],
           ",", 
          RowBox[{
          "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], 
         "]"}], "&"}], ",", "rz", ",", 
       RowBox[{"{", "2", "}"}]}], "]"}]}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"CellPrint", "[", 
     RowBox[{"Cell", "[", 
      RowBox[{
       RowBox[{
       "\"\<\:0420\:0430\:0441\:0441\:043c\:0430\:0442\:0440\:0438\:0432\:0430\
\:0435\:043c Re=\>\"", "<>", 
        RowBox[{"ToString", "[", "Rm", "]"}]}], ",", "\"\<Section\>\""}], 
      "]"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"len", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"Select", "[", 
       RowBox[{"len1", ",", 
        RowBox[{
         RowBox[{"0", "<", 
          RowBox[{
           RowBox[{"First", "[", "#", "]"}], "-", 
           RowBox[{"4", 
            RowBox[{"(", 
             RowBox[{"Rm", "-", "70"}], ")"}]}]}], "<", "40"}], "&"}]}], 
       "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"ltm", "=", 
     RowBox[{"First", "[", "len", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "mp", ",", "mdp", ",", "mvx", ",", "mdvx", ",", "mvy", ",", "mdvy", ",",
        "mvz", ",", "mdvz", ",", "qp", ",", "qvx", ",", "qvy", ",", "qvz"}], 
      "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Union", "[", 
        RowBox[{"Transpose", "[", 
         RowBox[{"{", 
          RowBox[{"ltm", ",", "#"}], "}"}], "]"}], "]"}], "&"}], "/@", 
      RowBox[{"Drop", "[", 
       RowBox[{"len", ",", "2"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"PrintRange", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
            "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "/", 
            "3"}], ",", 
           RowBox[{"#", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "&"}], "/@", 
        RowBox[{"(", 
         RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}]}], ",", "Automatic", 
       ",", 
       RowBox[{"PlotLabel", "\[Rule]", 
        SqrtBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}], 
          "\[LeftDoubleBracket]", 
          RowBox[{
           RowBox[{"-", "1"}], ",", 
           RowBox[{"-", "1"}]}], "\[RightDoubleBracket]"}]]}]}], "]"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"MapThread", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"PrintRange", "[", 
         RowBox[{"#1", ",", "Automatic", ",", 
          RowBox[{"PlotLabel", "\[Rule]", "#2"}]}], "]"}], "&"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
              "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", 
              RowBox[{"10", "^", 
               RowBox[{"Log", "[", 
                RowBox[{"10", ",", 
                 RowBox[{
                 "#", "\[LeftDoubleBracket]", "2", 
                  "\[RightDoubleBracket]"}]}], "]"}]}]}], "}"}], "&"}], ",", 
           RowBox[{"{", 
            RowBox[{"mvx", ",", "mvy", ",", "mvz", ",", "mp"}], "}"}], ",", 
           RowBox[{"{", "2", "}"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
          "\"\<Maximum of vr\>\"", ",", "\"\<Maximum of v\[Phi]\>\"", ",", 
           "\"\<Maximum of vz\>\"", ",", "\"\<Maximum of p\>\""}], "}"}]}], 
        "}"}]}], "]"}], "//", "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"$outp", "=", 
     RowBox[{"OpenWrite", "[", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"$Output", "=", 
     RowBox[{"{", "$outp", "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "maxvx", ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", 
       "mqvz"}], "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"ListLimit", "[", 
        RowBox[{
         RowBox[{"Take", "[", 
          RowBox[{
           RowBox[{"Smooth", "[", "#", "]"}], ",", 
           RowBox[{"-", 
            RowBox[{"Min", "[", 
             RowBox[{"500", ",", 
              RowBox[{"Length", "[", "mvy", "]"}]}], "]"}]}]}], "]"}], ",", 
         RowBox[{"CutValues", "\[Rule]", "False"}], ",", 
         RowBox[{"UseShift", "\[Rule]", "True"}], ",", 
         RowBox[{"ApproximationReport", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{"MaxError", ",", "HoelderError"}], "}"}]}]}], "]"}], "&"}],
       "/@", 
      RowBox[{"{", 
       RowBox[{
       "mvx", ",", "mvy", ",", "mvz", ",", "qvx", ",", "qvy", ",", "qvz"}], 
       "}"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Close", "[", "$outp", "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"$Output", "=", 
     RowBox[{"{", 
      RowBox[{"OutputStream", "[", 
       RowBox[{"\"\<stdout\>\"", ",", "1"}], "]"}], "}"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "dvx", ",", "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}], 
      "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", 
        RowBox[{"FromDigits", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Take", "[", 
            RowBox[{
             RowBox[{"First", "[", "#", "]"}], ",", "3"}], "]"}], ",", 
           RowBox[{"Last", "[", "#", "]"}]}], "}"}], "]"}], "]"}], "&"}], "/@", 
      RowBox[{"RealDigits", "[", 
       RowBox[{"First", "/@", 
        RowBox[{"ReadList", "[", 
         RowBox[{"First", "@", "$outp"}], "]"}]}], "]"}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"TableForm", "@", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "\"\<\>\"", ",", "\"\<vx\>\"", ",", "\"\<vy\>\"", ",", "\"\<vz\>\"", 
          ",", "\"\<Evx\>\"", ",", "\"\<Evy\>\"", ",", "\"\<Evz\>\""}], "}"}],
         ",", 
        RowBox[{"{", 
         RowBox[{
         "\"\<\:043f\:0440\:0435\:0434\:0435\:043b\:044b:\>\"", ",", "maxvx", 
          ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", "mqvz"}],
          "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "\"\<\:043e\:0448\:0438\:0431\:043a\:0438:\>\"", ",", "dvx", ",", 
          "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}], "}"}]}],
        "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"DeleteFile", "[", 
     RowBox[{"First", "@", "$outp"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"p", ",", "vx", ",", "vy", ",", "vz", ",", "nut"}], "}"}], "=", 
     RowBox[{"ReadTextSnapFile2D", "[", 
      RowBox[{
       RowBox[{"\"\<pumpflowRe\>\"", "<>", 
        RowBox[{"ToString", "[", "Rm", "]"}], "<>", "\"\<.snp\>\""}], ",", 
       "5", ",", 
       RowBox[{"PrintInfo", "\[Rule]", "Automatic"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<\:0421\:0434\:0432\:0438\:0433:\>\"", ",", 
      RowBox[{"shift", "=", 
       RowBox[{"Chop", "[", 
        RowBox[{"Block", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"fi", "=", 
             RowBox[{"Interpolation", "[", 
              RowBox[{"dif", "@", 
               RowBox[{"Transpose", "[", 
                RowBox[{"{", 
                 RowBox[{
                  RowBox[{
                  "coor", "\[LeftDoubleBracket]", "1", 
                   "\[RightDoubleBracket]"}], ",", 
                  RowBox[{"vy", "[", 
                   RowBox[{"[", 
                    RowBox[{"All", ",", 
                    RowBox[{"\[LeftFloor]", 
                    FractionBox[
                    RowBox[{"Length", "[", 
                    SuperscriptBox["vy", "\[Transpose]"], "]"}], "2"], 
                    "\[RightFloor]"}]}], "]"}], "]"}]}], "}"}], "]"}]}], 
              "]"}]}], ",", "t$"}], "}"}], ",", 
          RowBox[{"t$", "/.", 
           RowBox[{"FindRoot", "[", 
            RowBox[{
             RowBox[{"fi", "[", "t$", "]"}], ",", 
             RowBox[{"{", 
              RowBox[{"t$", ",", 
               RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}]}]}], "]"}], 
        "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"gr", "=", 
       RowBox[{"ListContourPlot", "[", 
        RowBox[{
         RowBox[{"GetInner", "[", 
          SuperscriptBox["vy", "\[Transpose]"], "]"}], ",", 
         RowBox[{"PlotLabel", "\[Rule]", 
          RowBox[{"\"\<\[Kappa]=\>\"", "<>", 
           RowBox[{"ToString", "[", "kappa", "]"}], "<>", "\"\<; Re=\>\"", "<>", 
           RowBox[{"ToString", "[", "Rm", "]"}]}]}], ",", 
         RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
         RowBox[{"Contours", "\[Rule]", "40"}]}], "]"}]}], ")"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Export", "[", 
     RowBox[{
      RowBox[{"\"\<pumpflowRe\>\"", "<>", 
       RowBox[{"ToString", "[", "Rm", "]"}], "<>", "\"\<v.gif\>\""}], ",", 
      "gr", ",", 
      RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"gr", "=", 
       RowBox[{"ListStreamPlot", "[", 
        RowBox[{
         RowBox[{"Transpose", "[", 
          RowBox[{
           RowBox[{"GetInner", "/@", 
            RowBox[{"{", 
             RowBox[{"vx", ",", "vz"}], "}"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"3", ",", "1", ",", "2"}], "}"}]}], "]"}], ",", 
         RowBox[{"StreamPoints", "\[Rule]", "Fine"}]}], "]"}]}], ")"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Export", "[", 
     RowBox[{
      RowBox[{"\"\<pumpflowRe\>\"", "<>", 
       RowBox[{"ToString", "[", "Rm", "]"}], "<>", "\"\<w.gif\>\""}], ",", 
      "gr", ",", 
      RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]z", "=", 
     RowBox[{
      RowBox[{"-", "M"}], " ", "vx"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Psi]r", "=", 
     RowBox[{"M", " ", "vz"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Psi]1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Last", "[", 
        RowBox[{"Transpose", "[", 
         RowBox[{"ListIntegration", "[", 
          RowBox[{"Transpose", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"Last", "[", "coor", "]"}], ",", "#"}], "}"}], "]"}], 
          "]"}], "]"}], "]"}], "&"}], "/@", "\[Psi]z"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]2", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Last", "[", 
         RowBox[{"Transpose", "[", 
          RowBox[{"ListIntegration", "[", 
           RowBox[{"Transpose", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{
              "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", "#"}], "}"}], "]"}], "]"}], "]"}], "]"}], "&"}], "/@", 
       RowBox[{"Transpose", "[", "\[Psi]r", "]"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Psi]1", "+", "\[Psi]2"}], ")"}], "/", "2"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
     "\"\<\:0444\:0443\:043d\:043a\:0446\:0438\:044f \:0442\:043e\:043a\:0430 \
\:0438 \:0435\:0451 \:043f\:043e\:0433\:0440\:0435\:0448\:043d\:043e\:0441\
\:0442\:044c:\>\"", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Max", "[", "#", "]"}], ",", 
          RowBox[{"Min", "[", "#", "]"}]}], "}"}], "&"}], "/@", 
       RowBox[{"{", 
        RowBox[{"\[Psi]", ",", 
         RowBox[{"\[Psi]2", "-", "\[Psi]1"}]}], "}"}]}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ListContourPlot", "[", 
      RowBox[{
       RowBox[{"Transpose", "[", "\[Psi]", "]"}], ",", 
       RowBox[{"Contours", "\[Rule]", "15"}], ",", 
       RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], "//", "Print"}], ";",
     "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
     "\"\<\:0417\:0430\:043f\:0438\:0441\:044b\:0432\:0430\:0435\:043c \:0440\
\:0435\:0437\:0443\:043b\:044c\:0442\:0430\:0442:\>\"", ",", 
      RowBox[{"putlist", "=", 
       RowBox[{
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"2011", ",", "3", ",", "2"}], "}"}], ",", 
          RowBox[{"{", "256", "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"Rm", ",", "kappa", ",", 
            RowBox[{"gamma", "=", 
             SuperscriptBox["10", 
              RowBox[{"-", "4"}]]}], ",", 
            RowBox[{"Fr", "=", "\"\<dp\>\""}], ",", 
            RowBox[{"Coriolis", "=", "0"}], ",", "shift", ",", "maxvx", ",", 
            "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", 
            "mqvx", ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", 
            "dqvz", ",", 
            RowBox[{"Max", "[", 
             RowBox[{"Abs", "[", "\[Psi]", "]"}], "]"}], ",", 
            RowBox[{"Max", "[", 
             RowBox[{"Abs", "[", 
              RowBox[{"\[Psi]2", "-", "\[Psi]1"}], "]"}], "]"}], ",", 
            "\"\<ErrScale 1e-4\>\""}], "}"}]}], "]"}], "/.", 
        RowBox[{"b_", "\[RuleDelayed]", 
         RowBox[{
          RowBox[{"CForm", "[", "b", "]"}], "/;", 
          RowBox[{"NumericQ", "[", "b", "]"}]}]}]}]}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"$resout", "=", 
     RowBox[{"OpenAppend", "[", 
      RowBox[{"$WD", "<>", "\"\<\\\\result.dat\>\""}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"SetOptions", "[", 
     RowBox[{"$resout", ",", 
      RowBox[{"PageWidth", "\[Rule]", "\[Infinity]"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"PutAppend", "[", 
     RowBox[{"putlist", ",", "$resout"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "gamma", ",", "Fr", ",", "Coriolis", ",", "shift", ",", "maxvx", ",", 
       "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", "mqvx",
        ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", "dqvz", ",", 
       "\[Psi]", ",", "\[Psi]2", ",", "\[Psi]1"}], "}"}], "=."}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Close", "[", "$resout", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"Rm", ",", "70", ",", "80", ",", "10"}], "}"}]}], 
  "\[IndentingNewLine]", "]"}]}], "Input",
 CellChangeTimes->{{3.424420539773*^9, 3.424420547672*^9}, {3.424420707809*^9,
    3.424420713389*^9}, {3.424420759052*^9, 3.4244207703389997`*^9}, 
   3.4244964698046*^9, {3.4244965404196*^9, 3.4244965907646*^9}, {
   3.424672039992*^9, 3.424672040644*^9}, {3.424685877091*^9, 
   3.424685898726*^9}, {3.4246860599700003`*^9, 3.424686089901*^9}, {
   3.4250226403256674`*^9, 3.4250226474017897`*^9}, {3.4250279418905106`*^9, 
   3.4250279748438053`*^9}, {3.425028351671485*^9, 3.4250283528966074`*^9}, {
   3.4420462538125*^9, 3.442046254796875*^9}, {3.442048608265625*^9, 
   3.44204860865625*^9}, {3.44205356553125*^9, 3.44205356628125*^9}, {
   3.4423126648125*^9, 3.442312680359375*^9}, {3.442318168484375*^9, 
   3.442318170390625*^9}, {3.4425015690810003`*^9, 3.4425015717539997`*^9}, 
   3.44250279999*^9, {3.442504533707*^9, 3.4425045570290003`*^9}, {
   3.442505197054*^9, 3.442505219318*^9}, {3.4425054209849997`*^9, 
   3.4425054225360003`*^9}, {3.442505601201*^9, 3.442505603276*^9}, {
   3.442505798258*^9, 3.442505800084*^9}, 3.4425065238059998`*^9, 
   3.4425065953199997`*^9, 3.442506774941*^9, {3.442507420257*^9, 
   3.442507436237*^9}, {3.4425080223050003`*^9, 3.442508042734*^9}, {
   3.442508112755*^9, 3.442508123516*^9}, {3.442508278344*^9, 
   3.442508279571*^9}, 3.442508681082*^9, {3.442508714898*^9, 
   3.4425087159779997`*^9}, {3.442509097354*^9, 3.442509115538*^9}, 
   3.4425091561*^9, 3.442509252835*^9, 3.442509310083*^9, 3.442509938058*^9, {
   3.442509980337*^9, 3.442509981443*^9}, {3.442510199034*^9, 
   3.442510247355*^9}, 3.4425110774*^9, {3.4425112942530003`*^9, 
   3.442511307491*^9}, {3.442511348616*^9, 3.442511349451*^9}, 
   3.4425115583310003`*^9, {3.4425116599820004`*^9, 3.442511660921*^9}, 
   3.442511717558*^9, {3.442512200615*^9, 3.4425122011619997`*^9}, {
   3.44251223701*^9, 3.442512238117*^9}, 3.442512555877*^9, 
   3.4425128875220003`*^9, 3.442514001566*^9, 3.44256471477*^9, 
   3.4426384774110003`*^9, 3.442639268377*^9, {3.442641432568*^9, 
   3.442641433236*^9}, {3.442641835356*^9, 3.4426418372720003`*^9}, 
   3.442642205058*^9, {3.442646594682*^9, 3.4426465954*^9}, {
   3.442647040994*^9, 3.4426470828789997`*^9}, {3.44264793402*^9, 
   3.442647934867*^9}, {3.4426488865690002`*^9, 3.442648886749*^9}, 
   3.442649348238*^9, {3.4426493890179996`*^9, 3.4426493899700003`*^9}, 
   3.44265027317*^9, {3.44265044346*^9, 3.4426504742539997`*^9}, {
   3.442650764574*^9, 3.4426507817720003`*^9}, {3.4428246783129997`*^9, 
   3.442824679712*^9}, {3.442825593979*^9, 3.442825611593*^9}, {
   3.442826304251*^9, 3.4428263061*^9}, 3.442826507415*^9, {
   3.4428307511610003`*^9, 3.442830752031*^9}, {3.44283107015*^9, 
   3.442831070665*^9}, {3.442831619543*^9, 3.442831620138*^9}, 
   3.4428346265030003`*^9, {3.4428351738570004`*^9, 3.4428351744309998`*^9}, {
   3.4428365579779997`*^9, 3.4428365583129997`*^9}, {3.442853403668*^9, 
   3.4428534040480003`*^9}, 3.44285369695*^9, {3.44584066375*^9, 
   3.445840664125*^9}, {3.4458407186875*^9, 3.4458407191875*^9}, {
   3.445840853984375*^9, 3.44584085521875*^9}, 3.445843044375*^9, {
   3.445843246828125*^9, 3.445843247234375*^9}, {3.445844257*^9, 
   3.445844257140625*^9}, {3.448352972538*^9, 3.4483529916005*^9}, {
   3.448694549109375*^9, 3.448694560546875*^9}, {3.448696574015625*^9, 
   3.448696575890625*^9}, {3.4486981549375*^9, 3.44869815709375*^9}, {
   3.448708845234375*^9, 3.44870884571875*^9}, {3.448709462359375*^9, 
   3.4487094626875*^9}, {3.44877703125*^9, 3.44877703184375*^9}, {
   3.4487774285625*^9, 3.44877742890625*^9}, {3.44878377809375*^9, 
   3.448783779171875*^9}, {3.448955335828125*^9, 3.448955336546875*^9}, {
   3.448955796*^9, 3.4489557978125*^9}, {3.448958334796875*^9, 
   3.448958337421875*^9}, {3.448958455515625*^9, 3.448958458015625*^9}, {
   3.44895854878125*^9, 3.448958552203125*^9}, 3.448958941453125*^9, {
   3.44895905065625*^9, 3.44895905071875*^9}, {3.4489591439375*^9, 
   3.448959146078125*^9}, {3.44895923565625*^9, 3.448959238578125*^9}, {
   3.448959301125*^9, 3.448959305109375*^9}, {3.448959353078125*^9, 
   3.448959355578125*^9}, {3.44895940759375*^9, 3.4489594100625*^9}, {
   3.448959975265625*^9, 3.448959979765625*^9}, {3.44896093871875*^9, 
   3.44896093884375*^9}, {3.44896115240625*^9, 3.448961156625*^9}, {
   3.44896353146875*^9, 3.44896353165625*^9}, {3.44896374259375*^9, 
   3.448963743015625*^9}, 3.449203645671875*^9, {3.4492118018125*^9, 
   3.4492118168125*^9}, {3.4492122759375*^9, 3.44921228315625*^9}, {
   3.44921680759375*^9, 3.449216843125*^9}, 3.449217611546875*^9, 
   3.44929218046875*^9, {3.449295434734375*^9, 3.44929546225*^9}, {
   3.449296950953125*^9, 3.4492969510625*^9}, {3.449297002046875*^9, 
   3.44929700490625*^9}, {3.449297185453125*^9, 3.449297187703125*^9}, {
   3.44929729690625*^9, 3.449297300375*^9}, {3.449297468234375*^9, 
   3.449297471046875*^9}, {3.449297553953125*^9, 3.4492975564375*^9}, {
   3.449297641921875*^9, 3.449297645671875*^9}, {3.449299009359375*^9, 
   3.449299014265625*^9}, {3.44982346515625*^9, 3.449823465421875*^9}, {
   3.449907878453125*^9, 3.449907888359375*^9}, {3.4499107680625*^9, 
   3.44991076996875*^9}, {3.4499959683125*^9, 3.4499959696875*^9}, {
   3.449996211296875*^9, 3.44999621390625*^9}, {3.449998746234375*^9, 
   3.449998750296875*^9}, {3.45007137253125*^9, 3.450071375453125*^9}, 
   3.45023583675*^9, 3.450236135609375*^9, {3.45023622740625*^9, 
   3.450236227453125*^9}, {3.45023656240625*^9, 3.4502365649375*^9}, {
   3.45041388878125*^9, 3.45041389415625*^9}, {3.45041413715625*^9, 
   3.450414137734375*^9}, {3.45041417090625*^9, 3.450414181234375*^9}, {
   3.450415025359375*^9, 3.4504150275625*^9}, {3.45041540059375*^9, 
   3.450415403328125*^9}, {3.450416030171875*^9, 3.45041603640625*^9}, {
   3.450416279515625*^9, 3.45041628184375*^9}, {3.450416388359375*^9, 
   3.450416390578125*^9}, {3.45042328396875*^9, 3.450423290859375*^9}, 
   3.458968236234375*^9, 3.4590535314375*^9, {3.459054455*^9, 
   3.45905450815625*^9}, {3.459061562640625*^9, 3.459061565078125*^9}, {
   3.4590616566875*^9, 3.459061658890625*^9}, {3.45906337284375*^9, 
   3.45906340934375*^9}, {3.4590637503125*^9, 3.459063758140625*^9}, {
   3.4590638070625*^9, 3.459063810453125*^9}, {3.45906389884375*^9, 
   3.45906390946875*^9}, {3.467008391978125*^9, 3.467008420759375*^9}, {
   3.467008476165625*^9, 3.467008476478125*^9}, {3.501441003135125*^9, 
   3.501441004916375*^9}, 3.5014827584375*^9, {3.501490815078125*^9, 
   3.50149081525*^9}, {3.50156001284375*^9, 3.501560013453125*^9}, {
   3.5015600833125*^9, 3.5015600836875*^9}, 3.501560358265625*^9, {
   3.505113046921875*^9, 3.50511306275*^9}, {3.5051180150625*^9, 
   3.50511801753125*^9}, {3.505118066703125*^9, 3.50511807196875*^9}, {
   3.505118118234375*^9, 3.505118136234375*^9}, 3.505118294375*^9, {
   3.50511879509375*^9, 3.505118796203125*^9}, {3.505118957703125*^9, 
   3.505118959125*^9}, {3.50511912059375*^9, 3.505119138953125*^9}, 
   3.505119262234375*^9, {3.50528246184375*^9, 3.50528247378125*^9}, {
   3.5052825175625*^9, 3.505282517671875*^9}, 3.505282779765625*^9, 
   3.505282938296875*^9, 3.50528306078125*^9, {3.505283458328125*^9, 
   3.50528348*^9}, {3.505291733640625*^9, 3.5052917338125*^9}, {
   3.50529179790625*^9, 3.50529180540625*^9}, 3.5052945143125*^9, 
   3.505295601703125*^9, 3.505295856453125*^9, 3.50529605653125*^9, 
   3.505296211390625*^9, {3.505296360078125*^9, 3.50529636015625*^9}, {
   3.505297056109375*^9, 3.505297056171875*^9}, 3.505297159703125*^9, 
   3.50529720196875*^9, 3.50529725765625*^9, {3.5052994333125*^9, 
   3.50529943346875*^9}, {3.505299467484375*^9, 3.505299512*^9}, 
   3.505299577765625*^9, 3.505299613703125*^9, 3.50529966634375*^9, 
   3.505299730640625*^9, {3.50535511421875*^9, 3.50535511865625*^9}, 
   3.50535541440625*^9, {3.505355467078125*^9, 3.505355477578125*^9}, {
   3.505357024921875*^9, 3.505357044921875*^9}, {3.5053570835*^9, 
   3.505357219109375*^9}, {3.505357338765625*^9, 3.5053575203125*^9}, 
   3.505357602359375*^9, {3.50535810834375*^9, 3.50535813946875*^9}, {
   3.505358394796875*^9, 3.5053584214375*^9}, {3.505358804375*^9, 
   3.50535881421875*^9}, {3.505358952203125*^9, 3.505358973203125*^9}, {
   3.505359078140625*^9, 3.505359122578125*^9}, {3.505359164984375*^9, 
   3.505359165109375*^9}, {3.505359248328125*^9, 3.505359249859375*^9}, {
   3.505359452953125*^9, 3.505359458109375*^9}, {3.505359490375*^9, 
   3.505359526125*^9}, {3.505400702890625*^9, 3.505400704640625*^9}, {
   3.505406893046875*^9, 3.5054068966875*^9}, {3.505442598890625*^9, 
   3.50544260890625*^9}, {3.5054451786875*^9, 3.505445186890625*^9}, {
   3.50544936884375*^9, 3.50544938628125*^9}, {3.5054497074375*^9, 
   3.505449731984375*^9}, {3.505449792859375*^9, 3.50544984690625*^9}, {
   3.505449993375*^9, 3.50545000253125*^9}, {3.505450057203125*^9, 
   3.50545006328125*^9}, {3.505450154*^9, 3.505450159546875*^9}, {
   3.50545035234375*^9, 3.505450357078125*^9}, {3.505450426*^9, 
   3.505450426859375*^9}, {3.50545066671875*^9, 3.505450670546875*^9}, {
   3.505450710234375*^9, 3.505450711984375*^9}, {3.505450840015625*^9, 
   3.505450840671875*^9}, {3.50545765009375*^9, 3.50545769615625*^9}, {
   3.505457950609375*^9, 3.505457950859375*^9}, {3.505458372125*^9, 
   3.5054583735*^9}, {3.505458766984375*^9, 3.505458771140625*^9}, {
   3.5054605455*^9, 3.5054605469375*^9}, {3.50546065390625*^9, 
   3.505460657*^9}, {3.50546209325*^9, 3.50546209721875*^9}, {
   3.505462560859375*^9, 3.5054625669375*^9}, {3.5054634454375*^9, 
   3.505463473203125*^9}, 3.505463539421875*^9, {3.505465609046875*^9, 
   3.505465621953125*^9}, {3.50546604915625*^9, 3.505466049203125*^9}, {
   3.505466238796875*^9, 3.50546624878125*^9}, {3.505468282171875*^9, 
   3.505468284140625*^9}, {3.505468323734375*^9, 3.505468332953125*^9}, {
   3.50546837353125*^9, 3.505468373609375*^9}, {3.50614766903275*^9, 
   3.50614767034525*^9}, 3.50614771290775*^9, {3.506148037392125*^9, 
   3.5061480480015*^9}, {3.5061484353765*^9, 3.506148446954625*^9}, {
   3.506149208798375*^9, 3.50614921097025*^9}, {3.506150064985875*^9, 
   3.506150069735875*^9}, {3.50632120590775*^9, 3.506321207439*^9}, {
   3.506664842142125*^9, 3.5066648721265*^9}, {3.5066649803765*^9, 
   3.50666498140775*^9}, {3.506665037314*^9, 3.50666513084525*^9}, {
   3.5069207033805*^9, 3.506920717271125*^9}, {3.506920807083625*^9, 
   3.50692081241175*^9}, {3.50692457416175*^9, 3.506924576396125*^9}, {
   3.50692468828675*^9, 3.506924699708625*^9}, {3.506924762708625*^9, 
   3.506924845271125*^9}, {3.50692504816175*^9, 3.5069250501305*^9}, {
   3.50726579209925*^9, 3.50726579341175*^9}, {3.50726597703675*^9, 
   3.507265978052375*^9}, {3.507266184583625*^9, 3.507266209677375*^9}, 
   3.507266330568*^9, 3.50726636109925*^9, {3.507266601958625*^9, 
   3.507266603239875*^9}, {3.50726668303675*^9, 3.507266684396125*^9}, {
   3.50726675103675*^9, 3.50726676891175*^9}, {3.507267430958625*^9, 
   3.507267431083625*^9}, {3.5072681642555*^9, 3.5072681852555*^9}, {
   3.5072687438805*^9, 3.50726875234925*^9}, {3.507269760271125*^9, 
   3.507269766396125*^9}, 3.507269807177375*^9, {3.5072706378805*^9, 
   3.507270644568*^9}, {3.50728151903675*^9, 3.507281525083625*^9}, 
   3.5072816558805*^9, {3.50728224541175*^9, 3.507282252364875*^9}, 
   3.50728237853675*^9, {3.50728290734925*^9, 3.507282908568*^9}, {
   3.5376861396983747`*^9, 3.5376861476202497`*^9}, {3.5376863665889997`*^9, 
   3.5376863674952497`*^9}, {3.5376864010577497`*^9, 
   3.5376864026358747`*^9}, {3.5376864660577497`*^9, 
   3.5376864758702497`*^9}, {3.5376865284327497`*^9, 
   3.5376865629483747`*^9}, {3.5376866103389997`*^9, 
   3.5376866107139997`*^9}, {3.53812934890625*^9, 3.538129365015625*^9}, {
   3.538129418140625*^9, 3.538129419546875*^9}, {3.53812949071875*^9, 
   3.53812949515625*^9}, {3.6242595882382812`*^9, 3.6242595885371094`*^9}, 
   3.624259621352539*^9, {3.626605842370117*^9, 3.626605843279297*^9}, {
   3.6266059091601562`*^9, 3.626605950883789*^9}, {3.6266092654453125`*^9, 
   3.6266092663203125`*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["\:041d\:0435\:0441\:043a\:043e\:043b\:044c\:043a\:043e \[Kappa] \:043f\
\:0440\:0438 \:043e\:0434\:043d\:043e\:043c Re", "Title",
 CellChangeTimes->{{3.5376859382139997`*^9, 3.5376859608858747`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"CellPrint", "[", 
   RowBox[{"Cell", "[", 
    RowBox[{
    "\"\<\:041c\:043e\:0436\:043d\:043e \:0441\:0442\:0438\:0440\:0430\:0442\
\:044c :) \>\"", ",", "\"\<Title\>\""}], "]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\:0420\:0430\:0441\:0441\:043c\:0430\:0442\:0440\:0438\:0432\:0430\
\:0435\:0442\:0441\:044f Re=\>\"", ",", 
    RowBox[{"Rm", "=", "200"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ReadTextSnapFile2D", "[", 
    RowBox[{"\"\<pumpflow(0.35).snp\>\"", ",", "5", ",", 
     RowBox[{"PrintInfo", "\[Rule]", "None"}]}], "]"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"for", " ", "ghost"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"coor", ",", "bounds", ",", "node"}], "}"}], "=", 
   RowBox[{"ReadAuxFiles2D", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"\"\<coord\>\"", ",", "\"\<node\>\""}], "}"}], ",", 
     RowBox[{"ShiftVector", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}]}], ",", 
     RowBox[{"ProcDistrib", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"4", ",", "4"}], "}"}]}], ",", 
     RowBox[{"NumberOfRefs", "\[Rule]", "2"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"kappastr", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"PadRight", "[", 
        RowBox[{
         RowBox[{"Characters", "[", "#", "]"}], ",", 
         RowBox[{"Max", "[", 
          RowBox[{
           RowBox[{"StringLength", "[", "#", "]"}], ",", "4"}], "]"}], ",", 
         "\"\<0\>\""}], "]"}], "&"}], "[", 
      RowBox[{"ToString", "[", "kappa", "]"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"CellPrint", "[", 
     RowBox[{"Cell", "[", 
      RowBox[{
       RowBox[{
       "\"\<\:0420\:0430\:0441\:0441\:043c\:0430\:0442\:0440\:0438\:0432\:0430\
\:0435\:043c \[Kappa]=\>\"", "<>", 
        RowBox[{"ToString", "[", "kappa", "]"}]}], ",", "\"\<Section\>\""}], 
      "]"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"M", "=", 
     RowBox[{"1", "+", 
      RowBox[{"kappa", " ", 
       RowBox[{
       "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}]}]}], 
    ";", "\[IndentingNewLine]", 
    RowBox[{"len1", "=", 
     RowBox[{"Import", "[", 
      RowBox[{
       RowBox[{"\"\<energy(\>\"", "<>", "kappastr", "<>", "\"\<).dat\>\""}], 
       ",", "\"\<Table\>\""}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"len", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"Select", "[", 
       RowBox[{"len1", ",", 
        RowBox[{
         RowBox[{"0", "<", 
          RowBox[{"First", "[", "#", "]"}]}], "&"}]}], "]"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"ltm", "=", 
     RowBox[{"First", "[", "len", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "mp", ",", "mdp", ",", "mvx", ",", "mdvx", ",", "mvy", ",", "mdvy", ",",
        "mvz", ",", "mdvz", ",", "qp", ",", "qvx", ",", "qvy", ",", "qvz"}], 
      "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Union", "[", 
        RowBox[{"Transpose", "[", 
         RowBox[{"{", 
          RowBox[{"ltm", ",", "#"}], "}"}], "]"}], "]"}], "&"}], "/@", 
      RowBox[{"Drop", "[", 
       RowBox[{"len", ",", "2"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"PrintRange", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
            "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "/", 
            "3"}], ",", 
           RowBox[{"#", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "&"}], "/@", 
        RowBox[{"(", 
         RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}]}], ",", "Automatic", 
       ",", 
       RowBox[{"PlotLabel", "\[Rule]", 
        SqrtBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}], 
          "\[LeftDoubleBracket]", 
          RowBox[{
           RowBox[{"-", "1"}], ",", 
           RowBox[{"-", "1"}]}], "\[RightDoubleBracket]"}]]}]}], "]"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"MapThread", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"PrintRange", "[", 
         RowBox[{"#1", ",", "Automatic", ",", 
          RowBox[{"PlotLabel", "\[Rule]", "#2"}]}], "]"}], "&"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
              "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", 
              RowBox[{"10", "^", 
               RowBox[{"Log", "[", 
                RowBox[{"10", ",", 
                 RowBox[{
                 "#", "\[LeftDoubleBracket]", "2", 
                  "\[RightDoubleBracket]"}]}], "]"}]}]}], "}"}], "&"}], ",", 
           RowBox[{"{", 
            RowBox[{"mvx", ",", "mvy", ",", "mvz", ",", "mp"}], "}"}], ",", 
           RowBox[{"{", "2", "}"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
          "\"\<Maximum of vr\>\"", ",", "\"\<Maximum of v\[Phi]\>\"", ",", 
           "\"\<Maximum of vz\>\"", ",", "\"\<Maximum of p\>\""}], "}"}]}], 
        "}"}]}], "]"}], "//", "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"$outp", "=", 
     RowBox[{"OpenWrite", "[", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"$Output", "=", 
     RowBox[{"{", "$outp", "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "maxvx", ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", 
       "mqvz"}], "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"ListLimit", "[", 
        RowBox[{
         RowBox[{"Take", "[", 
          RowBox[{
           RowBox[{"Smooth", "[", "#", "]"}], ",", 
           RowBox[{"-", 
            RowBox[{"Min", "[", 
             RowBox[{"500", ",", 
              RowBox[{
               RowBox[{"Length", "[", "mvy", "]"}], "-", "10"}]}], "]"}]}]}], 
          "]"}], ",", 
         RowBox[{"CutValues", "\[Rule]", "False"}], ",", 
         RowBox[{"UseShift", "\[Rule]", "True"}], ",", 
         RowBox[{"ApproximationReport", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{"MaxError", ",", "HoelderError"}], "}"}]}]}], "]"}], "&"}],
       "/@", 
      RowBox[{"{", 
       RowBox[{
       "mvx", ",", "mvy", ",", "mvz", ",", "qvx", ",", "qvy", ",", "qvz"}], 
       "}"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Close", "[", "$outp", "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"$Output", "=", 
     RowBox[{"{", 
      RowBox[{"OutputStream", "[", 
       RowBox[{"\"\<stdout\>\"", ",", "1"}], "]"}], "}"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "dvx", ",", "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}], 
      "}"}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"N", "[", 
        RowBox[{"FromDigits", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Take", "[", 
            RowBox[{
             RowBox[{"First", "[", "#", "]"}], ",", 
             RowBox[{"Min", "[", 
              RowBox[{"3", ",", 
               RowBox[{"Length", "[", 
                RowBox[{"First", "[", "#", "]"}], "]"}]}], "]"}]}], "]"}], 
           ",", 
           RowBox[{"Last", "[", "#", "]"}]}], "}"}], "]"}], "]"}], "&"}], "/@", 
      RowBox[{"RealDigits", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"First", "[", "#", "]"}], "&"}], "/@", 
        RowBox[{"ReadList", "[", 
         RowBox[{"First", "@", "$outp"}], "]"}]}], "]"}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"TableForm", "@", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "\"\<\>\"", ",", "\"\<vx\>\"", ",", "\"\<vy\>\"", ",", "\"\<vz\>\"", 
          ",", "\"\<Evx\>\"", ",", "\"\<Evy\>\"", ",", "\"\<Evz\>\""}], "}"}],
         ",", 
        RowBox[{"{", 
         RowBox[{
         "\"\<\:043f\:0440\:0435\:0434\:0435\:043b\:044b:\>\"", ",", "maxvx", 
          ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", "mqvz"}],
          "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "\"\<\:043e\:0448\:0438\:0431\:043a\:0438:\>\"", ",", "dvx", ",", 
          "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}], "}"}]}],
        "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"DeleteFile", "[", 
     RowBox[{"First", "@", "$outp"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"p", ",", "vx", ",", "vy", ",", "vz", ",", "nut"}], "}"}], "=", 
     RowBox[{"ReadTextSnapFile2D", "[", 
      RowBox[{
       RowBox[{"\"\<pumpflow(\>\"", "<>", "kappastr", "<>", "\"\<).snp\>\""}],
        ",", "5", ",", 
       RowBox[{"PrintInfo", "\[Rule]", "Automatic"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<\:0421\:0434\:0432\:0438\:0433:\>\"", ",", 
      RowBox[{"shift", "=", 
       RowBox[{"Chop", "[", 
        RowBox[{"Block", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"fi", "=", 
             RowBox[{"Interpolation", "[", 
              RowBox[{"dif", "@", 
               RowBox[{"Transpose", "[", 
                RowBox[{"{", 
                 RowBox[{
                  RowBox[{
                  "coor", "\[LeftDoubleBracket]", "1", 
                   "\[RightDoubleBracket]"}], ",", 
                  RowBox[{"vy", "[", 
                   RowBox[{"[", 
                    RowBox[{"All", ",", 
                    RowBox[{"\[LeftFloor]", 
                    FractionBox[
                    RowBox[{"Length", "[", 
                    SuperscriptBox["vy", "\[Transpose]"], "]"}], "2"], 
                    "\[RightFloor]"}]}], "]"}], "]"}]}], "}"}], "]"}]}], 
              "]"}]}], ",", "t$"}], "}"}], ",", 
          RowBox[{"t$", "/.", 
           RowBox[{"FindRoot", "[", 
            RowBox[{
             RowBox[{"fi", "[", "t$", "]"}], ",", 
             RowBox[{"{", 
              RowBox[{"t$", ",", 
               RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}]}]}], "]"}], 
        "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"gr", "=", 
       RowBox[{"ListContourPlot", "[", 
        RowBox[{
         RowBox[{"GetInner", "[", 
          SuperscriptBox["vy", "\[Transpose]"], "]"}], ",", 
         RowBox[{"PlotLabel", "\[Rule]", 
          RowBox[{"\"\<\[Kappa]=\>\"", "<>", 
           RowBox[{"ToString", "[", "kappa", "]"}], "<>", "\"\<; Re=\>\"", "<>", 
           RowBox[{"ToString", "[", "Rm", "]"}]}]}], ",", 
         RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
         RowBox[{"Contours", "\[Rule]", "40"}]}], "]"}]}], ")"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Export", "[", 
     RowBox[{
      RowBox[{"\"\<pumpflow\>\"", "<>", "kappastr", "<>", "\"\<v.gif\>\""}], 
      ",", "gr", ",", 
      RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"gr", "=", 
       RowBox[{"ListStreamPlot", "[", 
        RowBox[{
         RowBox[{"Transpose", "[", 
          RowBox[{
           RowBox[{"GetInner", "/@", 
            RowBox[{"{", 
             RowBox[{"vx", ",", "vz"}], "}"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"3", ",", "1", ",", "2"}], "}"}]}], "]"}], ",", 
         RowBox[{"StreamPoints", "\[Rule]", "Fine"}]}], "]"}]}], ")"}], "//", 
     "Print"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Export", "[", 
     RowBox[{
      RowBox[{"\"\<pumpflow\>\"", "<>", "kappastr", "<>", "\"\<w.gif\>\""}], 
      ",", "gr", ",", 
      RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]z", "=", 
     RowBox[{
      RowBox[{"-", "M"}], " ", "vx"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Psi]r", "=", 
     RowBox[{"M", " ", "vz"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"\[Psi]1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Last", "[", 
        RowBox[{"Transpose", "[", 
         RowBox[{"ListIntegration", "[", 
          RowBox[{"Transpose", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"Last", "[", "coor", "]"}], ",", "#"}], "}"}], "]"}], 
          "]"}], "]"}], "]"}], "&"}], "/@", "\[Psi]z"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]2", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Last", "[", 
         RowBox[{"Transpose", "[", 
          RowBox[{"ListIntegration", "[", 
           RowBox[{"Transpose", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{
              "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", "#"}], "}"}], "]"}], "]"}], "]"}], "]"}], "&"}], "/@", 
       RowBox[{"Transpose", "[", "\[Psi]r", "]"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"\[Psi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Psi]1", "+", "\[Psi]2"}], ")"}], "/", "2"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
     "\"\<\:0444\:0443\:043d\:043a\:0446\:0438\:044f \:0442\:043e\:043a\:0430 \
\:0438 \:0435\:0451 \:043f\:043e\:0433\:0440\:0435\:0448\:043d\:043e\:0441\
\:0442\:044c:\>\"", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Max", "[", "#", "]"}], ",", 
          RowBox[{"Min", "[", "#", "]"}]}], "}"}], "&"}], "/@", 
       RowBox[{"{", 
        RowBox[{"\[Psi]", ",", 
         RowBox[{"\[Psi]2", "-", "\[Psi]1"}]}], "}"}]}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ListContourPlot", "[", 
      RowBox[{
       RowBox[{"Transpose", "[", "\[Psi]", "]"}], ",", 
       RowBox[{"Contours", "\[Rule]", "15"}], ",", 
       RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], "//", "Print"}], ";",
     "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
     "\"\<\:0417\:0430\:043f\:0438\:0441\:044b\:0432\:0430\:0435\:043c \:0440\
\:0435\:0437\:0443\:043b\:044c\:0442\:0430\:0442:\>\"", ",", 
      RowBox[{"putlist", "=", 
       RowBox[{
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"2014", ",", "8", ",", "23"}], "}"}], ",", 
          RowBox[{"{", "256", "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"Rm", ",", "kappa", ",", 
            RowBox[{"gamma", "=", 
             SuperscriptBox["10", 
              RowBox[{"-", "4"}]]}], ",", 
            RowBox[{"Fr", "=", "\"\<1\>\""}], ",", 
            RowBox[{"Coriolis", "=", "0"}], ",", "shift", ",", "maxvx", ",", 
            "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", 
            "mqvx", ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", 
            "dqvz", ",", 
            RowBox[{"Max", "[", 
             RowBox[{"Abs", "[", "\[Psi]", "]"}], "]"}], ",", 
            RowBox[{"Max", "[", 
             RowBox[{"Abs", "[", 
              RowBox[{"\[Psi]2", "-", "\[Psi]1"}], "]"}], "]"}], ",", 
            "\"\<ErrScale 1e-4\>\""}], "}"}]}], "]"}], "/.", 
        RowBox[{"b_", "\[RuleDelayed]", 
         RowBox[{
          RowBox[{"CForm", "[", "b", "]"}], "/;", 
          RowBox[{"NumericQ", "[", "b", "]"}]}]}]}]}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"$resout", "=", 
     RowBox[{"OpenAppend", "[", 
      RowBox[{"$WD", "<>", "\"\<\\\\result.dat\>\""}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"SetOptions", "[", 
     RowBox[{"$resout", ",", 
      RowBox[{"PageWidth", "\[Rule]", "\[Infinity]"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"PutAppend", "[", 
     RowBox[{"putlist", ",", "$resout"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "gamma", ",", "Fr", ",", "Coriolis", ",", "shift", ",", "maxvx", ",", 
       "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", "mqvx",
        ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", "dqvz", ",", 
       "\[Psi]", ",", "\[Psi]2", ",", "\[Psi]1"}], "}"}], "=."}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Close", "[", "$resout", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"kappa", ",", "0.35", ",", "0.45", ",", "0.1"}], "}"}]}], 
  "\[IndentingNewLine]", "]"}]}], "Input",
 CellChangeTimes->{{3.424420539773*^9, 3.424420547672*^9}, {3.424420707809*^9,
    3.424420713389*^9}, {3.424420759052*^9, 3.4244207703389997`*^9}, 
   3.4244964698046*^9, {3.4244965404196*^9, 3.4244965907646*^9}, {
   3.424672039992*^9, 3.424672040644*^9}, {3.424685877091*^9, 
   3.424685898726*^9}, {3.4246860599700003`*^9, 3.424686089901*^9}, {
   3.4250226403256674`*^9, 3.4250226474017897`*^9}, {3.4250279418905106`*^9, 
   3.4250279748438053`*^9}, {3.425028351671485*^9, 3.4250283528966074`*^9}, {
   3.4420462538125*^9, 3.442046254796875*^9}, {3.442048608265625*^9, 
   3.44204860865625*^9}, {3.44205356553125*^9, 3.44205356628125*^9}, {
   3.4423126648125*^9, 3.442312680359375*^9}, {3.442318168484375*^9, 
   3.442318170390625*^9}, {3.4425015690810003`*^9, 3.4425015717539997`*^9}, 
   3.44250279999*^9, {3.442504533707*^9, 3.4425045570290003`*^9}, {
   3.442505197054*^9, 3.442505219318*^9}, {3.4425054209849997`*^9, 
   3.4425054225360003`*^9}, {3.442505601201*^9, 3.442505603276*^9}, {
   3.442505798258*^9, 3.442505800084*^9}, 3.4425065238059998`*^9, 
   3.4425065953199997`*^9, 3.442506774941*^9, {3.442507420257*^9, 
   3.442507436237*^9}, {3.4425080223050003`*^9, 3.442508042734*^9}, {
   3.442508112755*^9, 3.442508123516*^9}, {3.442508278344*^9, 
   3.442508279571*^9}, 3.442508681082*^9, {3.442508714898*^9, 
   3.4425087159779997`*^9}, {3.442509097354*^9, 3.442509115538*^9}, 
   3.4425091561*^9, 3.442509252835*^9, 3.442509310083*^9, 3.442509938058*^9, {
   3.442509980337*^9, 3.442509981443*^9}, {3.442510199034*^9, 
   3.442510247355*^9}, 3.4425110774*^9, {3.4425112942530003`*^9, 
   3.442511307491*^9}, {3.442511348616*^9, 3.442511349451*^9}, 
   3.4425115583310003`*^9, {3.4425116599820004`*^9, 3.442511660921*^9}, 
   3.442511717558*^9, {3.442512200615*^9, 3.4425122011619997`*^9}, {
   3.44251223701*^9, 3.442512238117*^9}, 3.442512555877*^9, 
   3.4425128875220003`*^9, 3.442514001566*^9, 3.44256471477*^9, 
   3.4426384774110003`*^9, 3.442639268377*^9, {3.442641432568*^9, 
   3.442641433236*^9}, {3.442641835356*^9, 3.4426418372720003`*^9}, 
   3.442642205058*^9, {3.442646594682*^9, 3.4426465954*^9}, {
   3.442647040994*^9, 3.4426470828789997`*^9}, {3.44264793402*^9, 
   3.442647934867*^9}, {3.4426488865690002`*^9, 3.442648886749*^9}, 
   3.442649348238*^9, {3.4426493890179996`*^9, 3.4426493899700003`*^9}, 
   3.44265027317*^9, {3.44265044346*^9, 3.4426504742539997`*^9}, {
   3.442650764574*^9, 3.4426507817720003`*^9}, {3.4428246783129997`*^9, 
   3.442824679712*^9}, {3.442825593979*^9, 3.442825611593*^9}, {
   3.442826304251*^9, 3.4428263061*^9}, 3.442826507415*^9, {
   3.4428307511610003`*^9, 3.442830752031*^9}, {3.44283107015*^9, 
   3.442831070665*^9}, {3.442831619543*^9, 3.442831620138*^9}, 
   3.4428346265030003`*^9, {3.4428351738570004`*^9, 3.4428351744309998`*^9}, {
   3.4428365579779997`*^9, 3.4428365583129997`*^9}, {3.442853403668*^9, 
   3.4428534040480003`*^9}, 3.44285369695*^9, {3.44584066375*^9, 
   3.445840664125*^9}, {3.4458407186875*^9, 3.4458407191875*^9}, {
   3.445840853984375*^9, 3.44584085521875*^9}, 3.445843044375*^9, {
   3.445843246828125*^9, 3.445843247234375*^9}, {3.445844257*^9, 
   3.445844257140625*^9}, {3.448352972538*^9, 3.4483529916005*^9}, {
   3.448694549109375*^9, 3.448694560546875*^9}, {3.448696574015625*^9, 
   3.448696575890625*^9}, {3.4486981549375*^9, 3.44869815709375*^9}, {
   3.448708845234375*^9, 3.44870884571875*^9}, {3.448709462359375*^9, 
   3.4487094626875*^9}, {3.44877703125*^9, 3.44877703184375*^9}, {
   3.4487774285625*^9, 3.44877742890625*^9}, {3.44878377809375*^9, 
   3.448783779171875*^9}, {3.448955335828125*^9, 3.448955336546875*^9}, {
   3.448955796*^9, 3.4489557978125*^9}, {3.448958334796875*^9, 
   3.448958337421875*^9}, {3.448958455515625*^9, 3.448958458015625*^9}, {
   3.44895854878125*^9, 3.448958552203125*^9}, 3.448958941453125*^9, {
   3.44895905065625*^9, 3.44895905071875*^9}, {3.4489591439375*^9, 
   3.448959146078125*^9}, {3.44895923565625*^9, 3.448959238578125*^9}, {
   3.448959301125*^9, 3.448959305109375*^9}, {3.448959353078125*^9, 
   3.448959355578125*^9}, {3.44895940759375*^9, 3.4489594100625*^9}, {
   3.448959975265625*^9, 3.448959979765625*^9}, {3.44896093871875*^9, 
   3.44896093884375*^9}, {3.44896115240625*^9, 3.448961156625*^9}, {
   3.44896353146875*^9, 3.44896353165625*^9}, {3.44896374259375*^9, 
   3.448963743015625*^9}, 3.449203645671875*^9, {3.4492118018125*^9, 
   3.4492118168125*^9}, {3.4492122759375*^9, 3.44921228315625*^9}, {
   3.44921680759375*^9, 3.449216843125*^9}, 3.449217611546875*^9, 
   3.44929218046875*^9, {3.449295434734375*^9, 3.44929546225*^9}, {
   3.449296950953125*^9, 3.4492969510625*^9}, {3.449297002046875*^9, 
   3.44929700490625*^9}, {3.449297185453125*^9, 3.449297187703125*^9}, {
   3.44929729690625*^9, 3.449297300375*^9}, {3.449297468234375*^9, 
   3.449297471046875*^9}, {3.449297553953125*^9, 3.4492975564375*^9}, {
   3.449297641921875*^9, 3.449297645671875*^9}, {3.449299009359375*^9, 
   3.449299014265625*^9}, {3.44982346515625*^9, 3.449823465421875*^9}, {
   3.449907878453125*^9, 3.449907888359375*^9}, {3.4499107680625*^9, 
   3.44991076996875*^9}, {3.4499959683125*^9, 3.4499959696875*^9}, {
   3.449996211296875*^9, 3.44999621390625*^9}, {3.449998746234375*^9, 
   3.449998750296875*^9}, {3.45007137253125*^9, 3.450071375453125*^9}, 
   3.45023583675*^9, 3.450236135609375*^9, {3.45023622740625*^9, 
   3.450236227453125*^9}, {3.45023656240625*^9, 3.4502365649375*^9}, {
   3.45041388878125*^9, 3.45041389415625*^9}, {3.45041413715625*^9, 
   3.450414137734375*^9}, {3.45041417090625*^9, 3.450414181234375*^9}, {
   3.450415025359375*^9, 3.4504150275625*^9}, {3.45041540059375*^9, 
   3.450415403328125*^9}, {3.450416030171875*^9, 3.45041603640625*^9}, {
   3.450416279515625*^9, 3.45041628184375*^9}, {3.450416388359375*^9, 
   3.450416390578125*^9}, {3.45042328396875*^9, 3.450423290859375*^9}, 
   3.458968236234375*^9, 3.4590535314375*^9, {3.459054455*^9, 
   3.45905450815625*^9}, {3.459061562640625*^9, 3.459061565078125*^9}, {
   3.4590616566875*^9, 3.459061658890625*^9}, {3.45906337284375*^9, 
   3.45906340934375*^9}, {3.4590637503125*^9, 3.459063758140625*^9}, {
   3.4590638070625*^9, 3.459063810453125*^9}, {3.45906389884375*^9, 
   3.45906390946875*^9}, {3.467008391978125*^9, 3.467008420759375*^9}, {
   3.467008476165625*^9, 3.467008476478125*^9}, {3.501441003135125*^9, 
   3.501441004916375*^9}, 3.5014827584375*^9, {3.501490815078125*^9, 
   3.50149081525*^9}, {3.50156001284375*^9, 3.501560013453125*^9}, {
   3.5015600833125*^9, 3.5015600836875*^9}, 3.501560358265625*^9, {
   3.505113046921875*^9, 3.50511306275*^9}, {3.5051180150625*^9, 
   3.50511801753125*^9}, {3.505118066703125*^9, 3.50511807196875*^9}, {
   3.505118118234375*^9, 3.505118136234375*^9}, 3.505118294375*^9, {
   3.50511879509375*^9, 3.505118796203125*^9}, {3.505118957703125*^9, 
   3.505118959125*^9}, {3.50511912059375*^9, 3.505119138953125*^9}, 
   3.505119262234375*^9, {3.50528246184375*^9, 3.50528247378125*^9}, {
   3.5052825175625*^9, 3.505282517671875*^9}, 3.505282779765625*^9, 
   3.505282938296875*^9, 3.50528306078125*^9, {3.505283458328125*^9, 
   3.50528348*^9}, {3.505291733640625*^9, 3.5052917338125*^9}, {
   3.50529179790625*^9, 3.50529180540625*^9}, 3.5052945143125*^9, 
   3.505295601703125*^9, 3.505295856453125*^9, 3.50529605653125*^9, 
   3.505296211390625*^9, {3.505296360078125*^9, 3.50529636015625*^9}, {
   3.505297056109375*^9, 3.505297056171875*^9}, 3.505297159703125*^9, 
   3.50529720196875*^9, 3.50529725765625*^9, {3.5052994333125*^9, 
   3.50529943346875*^9}, {3.505299467484375*^9, 3.505299512*^9}, 
   3.505299577765625*^9, 3.505299613703125*^9, 3.50529966634375*^9, 
   3.505299730640625*^9, {3.50535511421875*^9, 3.50535511865625*^9}, 
   3.50535541440625*^9, {3.505355467078125*^9, 3.505355477578125*^9}, {
   3.505357024921875*^9, 3.505357044921875*^9}, {3.5053570835*^9, 
   3.505357219109375*^9}, {3.505357338765625*^9, 3.5053575203125*^9}, 
   3.505357602359375*^9, {3.50535810834375*^9, 3.50535813946875*^9}, {
   3.505358394796875*^9, 3.5053584214375*^9}, {3.505358804375*^9, 
   3.50535881421875*^9}, {3.505358952203125*^9, 3.505358973203125*^9}, {
   3.505359078140625*^9, 3.505359122578125*^9}, {3.505359164984375*^9, 
   3.505359165109375*^9}, {3.505359248328125*^9, 3.505359249859375*^9}, {
   3.505359452953125*^9, 3.505359458109375*^9}, {3.505359490375*^9, 
   3.505359526125*^9}, {3.505400702890625*^9, 3.505400704640625*^9}, {
   3.505406893046875*^9, 3.5054068966875*^9}, {3.505442598890625*^9, 
   3.50544260890625*^9}, {3.5054451786875*^9, 3.505445186890625*^9}, {
   3.50544936884375*^9, 3.50544938628125*^9}, {3.5054497074375*^9, 
   3.505449731984375*^9}, {3.505449792859375*^9, 3.50544984690625*^9}, {
   3.505449993375*^9, 3.50545000253125*^9}, {3.505450057203125*^9, 
   3.50545006328125*^9}, {3.505450154*^9, 3.505450159546875*^9}, {
   3.50545035234375*^9, 3.505450357078125*^9}, {3.505450426*^9, 
   3.505450426859375*^9}, {3.50545066671875*^9, 3.505450670546875*^9}, {
   3.505450710234375*^9, 3.505450711984375*^9}, {3.505450840015625*^9, 
   3.505450840671875*^9}, {3.50545765009375*^9, 3.50545769615625*^9}, {
   3.505457950609375*^9, 3.505457950859375*^9}, {3.505458372125*^9, 
   3.5054583735*^9}, {3.505458766984375*^9, 3.505458771140625*^9}, {
   3.5054605455*^9, 3.5054605469375*^9}, {3.50546065390625*^9, 
   3.505460657*^9}, {3.50546209325*^9, 3.50546209721875*^9}, {
   3.505462560859375*^9, 3.5054625669375*^9}, {3.5054634454375*^9, 
   3.505463473203125*^9}, 3.505463539421875*^9, {3.505465609046875*^9, 
   3.505465621953125*^9}, {3.50546604915625*^9, 3.505466049203125*^9}, {
   3.505466238796875*^9, 3.50546624878125*^9}, {3.505468282171875*^9, 
   3.505468284140625*^9}, {3.505468323734375*^9, 3.505468332953125*^9}, {
   3.50546837353125*^9, 3.505468373609375*^9}, {3.50614766903275*^9, 
   3.50614767034525*^9}, 3.50614771290775*^9, {3.506148037392125*^9, 
   3.5061480480015*^9}, {3.5061484353765*^9, 3.506148446954625*^9}, {
   3.506149208798375*^9, 3.50614921097025*^9}, {3.506150064985875*^9, 
   3.506150069735875*^9}, {3.50632120590775*^9, 3.506321207439*^9}, {
   3.506664842142125*^9, 3.5066648721265*^9}, {3.5066649803765*^9, 
   3.50666498140775*^9}, {3.506665037314*^9, 3.50666513084525*^9}, {
   3.5069207033805*^9, 3.506920717271125*^9}, {3.506920807083625*^9, 
   3.50692081241175*^9}, {3.50692457416175*^9, 3.506924576396125*^9}, {
   3.50692468828675*^9, 3.506924699708625*^9}, {3.506924762708625*^9, 
   3.506924845271125*^9}, {3.50692504816175*^9, 3.5069250501305*^9}, {
   3.50726579209925*^9, 3.50726579341175*^9}, {3.50726597703675*^9, 
   3.507265978052375*^9}, {3.507266184583625*^9, 3.507266209677375*^9}, 
   3.507266330568*^9, 3.50726636109925*^9, {3.507266601958625*^9, 
   3.507266603239875*^9}, {3.50726668303675*^9, 3.507266684396125*^9}, {
   3.50726675103675*^9, 3.50726676891175*^9}, {3.507267430958625*^9, 
   3.507267431083625*^9}, {3.5072681642555*^9, 3.5072681852555*^9}, {
   3.5072687438805*^9, 3.50726875234925*^9}, {3.507269760271125*^9, 
   3.507269766396125*^9}, 3.507269807177375*^9, {3.5072706378805*^9, 
   3.507270644568*^9}, {3.50728151903675*^9, 3.507281525083625*^9}, 
   3.5072816558805*^9, {3.50728224541175*^9, 3.507282252364875*^9}, 
   3.50728237853675*^9, {3.50728290734925*^9, 3.507282908568*^9}, 
   3.536049312431506*^9, {3.536049357744006*^9, 3.536049433244006*^9}, {
   3.5360494901658816`*^9, 3.536049578556506*^9}, 3.536049772212756*^9, {
   3.5360499023846316`*^9, 3.5360499034158816`*^9}, 3.5360499380408816`*^9, {
   3.536050075337756*^9, 3.5360500911971316`*^9}, {3.538384632640625*^9, 
   3.5383846836875*^9}, {3.53838471803125*^9, 3.53838473909375*^9}, 
   3.538384857234375*^9, {3.538385025984375*^9, 3.53838503109375*^9}, {
   3.5383851689375*^9, 3.5383851760625*^9}, {3.538385271234375*^9, 
   3.538385333921875*^9}, 3.53838547653125*^9, 3.53838555525*^9, {
   3.538385846125*^9, 3.538385859828125*^9}, {3.538385894515625*^9, 
   3.538385896859375*^9}, {3.538385930515625*^9, 3.538385961*^9}, 
   3.5383859910625*^9, {3.53838602515625*^9, 3.538386045453125*^9}, {
   3.53838616103125*^9, 3.538386161296875*^9}, {3.53838621303125*^9, 
   3.5383862245*^9}, {3.538386368375*^9, 3.5383863688125*^9}, {
   3.53838657903125*^9, 3.538386595921875*^9}, {3.5387380713125*^9, 
   3.5387380926875*^9}, {3.53873829765625*^9, 3.53873829809375*^9}, {
   3.539334962817375*^9, 3.53933496405175*^9}, {3.5998130324873123`*^9, 
   3.599813065052742*^9}, {3.599813729537117*^9, 3.5998137296904373`*^9}, {
   3.5998160101503983`*^9, 3.599816018771492*^9}, {3.5998949365400467`*^9, 
   3.599894955188484*^9}, 3.6036179341128206`*^9, 3.604213893998047*^9, {
   3.604213926665039*^9, 3.604213931001953*^9}, {3.604301955510742*^9, 
   3.604301998678711*^9}, {3.6052606314778123`*^9, 3.6052606367278123`*^9}, {
   3.6052607123420706`*^9, 3.6052607126877737`*^9}, {3.6176898697851405`*^9, 
   3.617689925326156*^9}, {3.617690245195297*^9, 3.617690246589828*^9}, {
   3.6176910863535*^9, 3.617691087021469*^9}, 3.6234821000664062`*^9, {
   3.6234821307978516`*^9, 3.623482134913086*^9}, {3.623482186390625*^9, 
   3.6234821902373047`*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["(\[Kappa], Re) \:0438\:0437 \:0444\:0430\:0439\:043b\:0430", "Title",
 CellChangeTimes->{{3.5376859382139997`*^9, 3.5376859608858747`*^9}, {
  3.604213212140625*^9, 3.604213217602539*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"CellPrint", "[", 
   RowBox[{"Cell", "[", 
    RowBox[{
    "\"\<\:041c\:043e\:0436\:043d\:043e \:0441\:0442\:0438\:0440\:0430\:0442\
\:044c :) \>\"", ",", "\"\<Title\>\""}], "]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\:0412 \:0444\:0430\:0439\:043b\:0435 \:0431\:044b\:043b\:0438 \:0437\
\:0430\:0434\:0430\:043d\:044b \:043f\:0430\:0440\:0430\:043c\:0435\:0442\
\:0440\:044b: \>\"", ",", 
    RowBox[{"First", "@", 
     RowBox[{"Import", "[", "\"\<pumpflow.params\>\"", "]"}]}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pars", "=", 
   RowBox[{"Import", "[", "\"\<pumpflow.done\>\"", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<\:041d\:0430\:0439\:0434\:0435\:043d\:043e \>\"", ",", 
   RowBox[{"Length", "[", "pars", "]"}], ",", 
   "\"\< \:043d\:0430\:0431\:043e\:0440\:043e\:0432 \:043f\:0430\:0440\:0430\
\:043c\:0435\:0442\:0440\:043e\:0432:\>\"", ",", 
   RowBox[{"TableForm", "[", "pars", "]"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ReadTextSnapFile2D", "[", 
   RowBox[{"\"\<pumpflow(kappa=0.50,Re=80).snp\>\"", ",", "5", ",", 
    RowBox[{"PrintInfo", "\[Rule]", "None"}]}], "]"}], ";", 
  RowBox[{"(*", 
   RowBox[{"for", " ", "ghost"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{"kappa", ",", "Rm", ",", "Ttot"}], "}"}], "=", 
      RowBox[{"Drop", "[", 
       RowBox[{
        RowBox[{"pars", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"{", "3", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"kappastr", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"PadRight", "[", 
         RowBox[{
          RowBox[{"Characters", "[", "#", "]"}], ",", 
          RowBox[{"Max", "[", 
           RowBox[{
            RowBox[{"StringLength", "[", "#", "]"}], ",", "4"}], "]"}], ",", 
          "\"\<0\>\""}], "]"}], "&"}], "[", 
       RowBox[{"ToString", "[", "kappa", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"parstr", "=", 
      RowBox[{"\"\<kappa=\>\"", "<>", "kappastr", "<>", "\"\<,Re=\>\"", "<>", 
       RowBox[{"ToString", "[", "Rm", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"CellPrint", "[", 
      RowBox[{"Cell", "[", 
       RowBox[{
        RowBox[{
        "\"\<\:0420\:0430\:0441\:0441\:043c\:0430\:0442\:0440\:0438\:0432\
\:0430\:0435\:043c \[Kappa]=\>\"", "<>", 
         RowBox[{"ToString", "[", "kappa", "]"}], "<>", "\"\<, Re=\>\"", "<>", 
         RowBox[{"ToString", "[", "Rm", "]"}]}], ",", "\"\<Section\>\""}], 
       "]"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{
      "\"\<\:0421\:0442\:0440\:043e\:043a\:0430 \:0434\:043b\:044f \:0444\
\:0430\:0439\:043b\:043e\:0432: (\>\"", ",", "parstr", ",", "\"\<)\>\""}], 
      "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"coor", ",", "bounds", ",", "node"}], "}"}], "=", 
      RowBox[{"ReadAuxFiles2D", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<coord\>\"", ",", "\"\<node\>\""}], "}"}], ",", 
        RowBox[{"ShiftVector", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{"0", ",", "0"}], "}"}]}], ",", 
        RowBox[{"ProcDistrib", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{"4", ",", "4"}], "}"}]}], ",", 
        RowBox[{"NumberOfRefs", "\[Rule]", "2"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M", "=", 
      RowBox[{"1", "+", 
       RowBox[{"kappa", " ", 
        RowBox[{
        "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}]}]}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"len1", "=", 
      RowBox[{"Import", "[", 
       RowBox[{
        RowBox[{"\"\<energy(\>\"", "<>", "parstr", "<>", "\"\<).dat\>\""}], 
        ",", "\"\<Table\>\""}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"len", "=", 
      RowBox[{"Transpose", "[", 
       RowBox[{"Select", "[", 
        RowBox[{"len1", ",", 
         RowBox[{
          RowBox[{"0", "<", 
           RowBox[{"First", "[", "#", "]"}]}], "&"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"ltm", "=", 
      RowBox[{"First", "[", "len", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "mp", ",", "mdp", ",", "mvx", ",", "mdvx", ",", "mvy", ",", "mdvy", 
        ",", "mvz", ",", "mdvz", ",", "qp", ",", "qvx", ",", "qvy", ",", 
        "qvz"}], "}"}], "=", 
      RowBox[{
       RowBox[{
        RowBox[{"Union", "[", 
         RowBox[{"Transpose", "[", 
          RowBox[{"{", 
           RowBox[{"ltm", ",", "#"}], "}"}], "]"}], "]"}], "&"}], "/@", 
       RowBox[{"Drop", "[", 
        RowBox[{"len", ",", "2"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"PrintRange", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{
             "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "/",
              "3"}], ",", 
            RowBox[{"#", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "&"}], "/@", 
         RowBox[{"(", 
          RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}]}], ",", "Automatic", 
        ",", 
        RowBox[{"PlotLabel", "\[Rule]", 
         SqrtBox[
          RowBox[{
           RowBox[{"(", 
            RowBox[{"qvx", "+", "qvy", "+", "qvz"}], ")"}], 
           "\[LeftDoubleBracket]", 
           RowBox[{
            RowBox[{"-", "1"}], ",", 
            RowBox[{"-", "1"}]}], "\[RightDoubleBracket]"}]]}]}], "]"}], "//",
       "Print"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"MapThread", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"PrintRange", "[", 
          RowBox[{"#1", ",", "Automatic", ",", 
           RowBox[{"PlotLabel", "\[Rule]", "#2"}]}], "]"}], "&"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Map", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{
               "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
               ",", 
               RowBox[{"10", "^", 
                RowBox[{"Log", "[", 
                 RowBox[{"10", ",", 
                  RowBox[{
                  "#", "\[LeftDoubleBracket]", "2", 
                   "\[RightDoubleBracket]"}]}], "]"}]}]}], "}"}], "&"}], ",", 
            RowBox[{"{", 
             RowBox[{"mvx", ",", "mvy", ",", "mvz", ",", "mp"}], "}"}], ",", 
            RowBox[{"{", "2", "}"}]}], "]"}], ",", 
          RowBox[{"{", 
           RowBox[{
           "\"\<Maximum of vr\>\"", ",", "\"\<Maximum of v\[Phi]\>\"", ",", 
            "\"\<Maximum of vz\>\"", ",", "\"\<Maximum of p\>\""}], "}"}]}], 
         "}"}]}], "]"}], "//", "Print"}], ";", "\[IndentingNewLine]", 
     RowBox[{"$outp", "=", 
      RowBox[{"OpenWrite", "[", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"$Output", "=", 
      RowBox[{"{", "$outp", "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "maxvx", ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", 
        "mqvz"}], "}"}], "=", 
      RowBox[{
       RowBox[{
        RowBox[{"ListLimit", "[", 
         RowBox[{
          RowBox[{"Smooth", "[", 
           RowBox[{"Take", "[", 
            RowBox[{"#", ",", 
             RowBox[{"-", 
              RowBox[{"Min", "[", 
               RowBox[{"5000", ",", 
                RowBox[{
                 RowBox[{"Length", "[", "mvy", "]"}], "-", "10"}]}], 
               "]"}]}]}], "]"}], "]"}], ",", 
          RowBox[{"CutValues", "\[Rule]", "False"}], ",", 
          RowBox[{"UseShift", "\[Rule]", "True"}], ",", 
          RowBox[{"ApproximationReport", "\[Rule]", 
           RowBox[{"{", 
            RowBox[{"MaxError", ",", "HoelderError"}], "}"}]}]}], "]"}], 
        "&"}], "/@", 
       RowBox[{"{", 
        RowBox[{
        "mvx", ",", "mvy", ",", "mvz", ",", "qvx", ",", "qvy", ",", "qvz"}], 
        "}"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Close", "[", "$outp", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"$Output", "=", 
      RowBox[{"{", 
       RowBox[{"OutputStream", "[", 
        RowBox[{"\"\<stdout\>\"", ",", "1"}], "]"}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "dvx", ",", "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}],
        "}"}], "=", 
      RowBox[{
       RowBox[{
        RowBox[{"N", "[", 
         RowBox[{"FromDigits", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Take", "[", 
             RowBox[{
              RowBox[{"First", "[", "#", "]"}], ",", 
              RowBox[{"Min", "[", 
               RowBox[{"3", ",", 
                RowBox[{"Length", "[", 
                 RowBox[{"First", "[", "#", "]"}], "]"}]}], "]"}]}], "]"}], 
            ",", 
            RowBox[{"Last", "[", "#", "]"}]}], "}"}], "]"}], "]"}], "&"}], "/@", 
       RowBox[{"RealDigits", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"First", "[", "#", "]"}], "&"}], "/@", 
         RowBox[{"ReadList", "[", 
          RowBox[{"First", "@", "$outp"}], "]"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"TableForm", "@", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
          "\"\<\>\"", ",", "\"\<vx\>\"", ",", "\"\<vy\>\"", ",", "\"\<vz\>\"",
            ",", "\"\<Evx\>\"", ",", "\"\<Evy\>\"", ",", "\"\<Evz\>\""}], 
          "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
          "\"\<\:043f\:0440\:0435\:0434\:0435\:043b\:044b:\>\"", ",", "maxvx",
            ",", "maxvy", ",", "maxvz", ",", "mqvx", ",", "mqvy", ",", 
           "mqvz"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
          "\"\<\:043e\:0448\:0438\:0431\:043a\:0438:\>\"", ",", "dvx", ",", 
           "dvy", ",", "dvz", ",", "dqvx", ",", "dqvy", ",", "dqvz"}], 
          "}"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"DeleteFile", "[", 
      RowBox[{"First", "@", "$outp"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"p", ",", "vx", ",", "vy", ",", "vz", ",", "nut"}], "}"}], "=", 
      RowBox[{"ReadTextSnapFile2D", "[", 
       RowBox[{
        RowBox[{"\"\<pumpflow(\>\"", "<>", "parstr", "<>", "\"\<).snp\>\""}], 
        ",", "5", ",", 
        RowBox[{"PrintInfo", "\[Rule]", "Automatic"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"\"\<\:0421\:0434\:0432\:0438\:0433:\>\"", ",", 
       RowBox[{"shift", "=", 
        RowBox[{"Chop", "[", 
         RowBox[{"Block", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"fi", "=", 
              RowBox[{"Interpolation", "[", 
               RowBox[{"dif", "@", 
                RowBox[{"Transpose", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{
                   "coor", "\[LeftDoubleBracket]", "1", 
                    "\[RightDoubleBracket]"}], ",", 
                   RowBox[{"vy", "[", 
                    RowBox[{"[", 
                    RowBox[{"All", ",", 
                    RowBox[{"\[LeftFloor]", 
                    FractionBox[
                    RowBox[{"Length", "[", 
                    SuperscriptBox["vy", "\[Transpose]"], "]"}], "2"], 
                    "\[RightFloor]"}]}], "]"}], "]"}]}], "}"}], "]"}]}], 
               "]"}]}], ",", "t$"}], "}"}], ",", 
           RowBox[{"t$", "/.", 
            RowBox[{"FindRoot", "[", 
             RowBox[{
              RowBox[{"fi", "[", "t$", "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"t$", ",", 
                RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}]}]}], "]"}], 
         "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"gr", "=", 
        RowBox[{"ListContourPlot", "[", 
         RowBox[{
          RowBox[{"GetInner", "[", 
           SuperscriptBox["vy", "\[Transpose]"], "]"}], ",", 
          RowBox[{"PlotLabel", "\[Rule]", 
           RowBox[{"\"\<\[Kappa]=\>\"", "<>", 
            RowBox[{"ToString", "[", "kappa", "]"}], "<>", "\"\<; Re=\>\"", "<>", 
            RowBox[{"ToString", "[", "Rm", "]"}]}]}], ",", 
          RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
          RowBox[{"Contours", "\[Rule]", "40"}]}], "]"}]}], ")"}], "//", 
      "Print"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"\"\<pumpflow(\>\"", "<>", "parstr", "<>", "\"\<)v.gif\>\""}], 
       ",", "gr", ",", 
       RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
       RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"gr", "=", 
        RowBox[{"ListStreamPlot", "[", 
         RowBox[{
          RowBox[{"Transpose", "[", 
           RowBox[{
            RowBox[{"GetInner", "/@", 
             RowBox[{"{", 
              RowBox[{"vx", ",", "vz"}], "}"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"3", ",", "1", ",", "2"}], "}"}]}], "]"}], ",", 
          RowBox[{"StreamPoints", "\[Rule]", "Fine"}]}], "]"}]}], ")"}], "//",
       "Print"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"\"\<pumpflow(\>\"", "<>", "parstr", "<>", "\"\<)w.gif\>\""}], 
       ",", "gr", ",", 
       RowBox[{"ImageResolution", "\[Rule]", "80"}], ",", 
       RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"\[Psi]z", "=", 
      RowBox[{
       RowBox[{"-", "M"}], " ", "vx"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"\[Psi]r", "=", 
      RowBox[{"M", " ", "vz"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"\[Psi]1", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"Last", "[", 
         RowBox[{"Transpose", "[", 
          RowBox[{"ListIntegration", "[", 
           RowBox[{"Transpose", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Last", "[", "coor", "]"}], ",", "#"}], "}"}], "]"}], 
           "]"}], "]"}], "]"}], "&"}], "/@", "\[Psi]z"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"\[Psi]2", "=", 
      RowBox[{"Transpose", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Last", "[", 
          RowBox[{"Transpose", "[", 
           RowBox[{"ListIntegration", "[", 
            RowBox[{"Transpose", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{
               "coor", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}],
                ",", "#"}], "}"}], "]"}], "]"}], "]"}], "]"}], "&"}], "/@", 
        RowBox[{"Transpose", "[", "\[Psi]r", "]"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"\[Psi]", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Psi]1", "+", "\[Psi]2"}], ")"}], "/", "2"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{
      "\"\<\:0444\:0443\:043d\:043a\:0446\:0438\:044f \
\:0442\:043e\:043a\:0430 \:0438 \:0435\:0451 \:043f\:043e\:0433\:0440\:0435\
\:0448\:043d\:043e\:0441\:0442\:044c:\>\"", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Max", "[", "#", "]"}], ",", 
           RowBox[{"Min", "[", "#", "]"}]}], "}"}], "&"}], "/@", 
        RowBox[{"{", 
         RowBox[{"\[Psi]", ",", 
          RowBox[{"\[Psi]2", "-", "\[Psi]1"}]}], "}"}]}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"ListContourPlot", "[", 
       RowBox[{
        RowBox[{"Transpose", "[", "\[Psi]", "]"}], ",", 
        RowBox[{"Contours", "\[Rule]", "15"}], ",", 
        RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], "//", "Print"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{
      "\"\<\:0417\:0430\:043f\:0438\:0441\:044b\:0432\:0430\:0435\:043c \
\:0440\:0435\:0437\:0443\:043b\:044c\:0442\:0430\:0442:\>\"", ",", 
       RowBox[{"putlist", "=", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"2014", ",", "12", ",", "3"}], "}"}], ",", 
           RowBox[{"{", "256", "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"Rm", ",", "kappa", ",", 
             RowBox[{"gamma", "=", 
              SuperscriptBox["10", 
               RowBox[{"-", "4"}]]}], ",", 
             RowBox[{"Fr", "=", "\"\<r\>\""}], ",", 
             RowBox[{"Coriolis", "=", "0"}], ",", "shift", ",", "maxvx", ",", 
             "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", 
             "mqvx", ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", 
             "dqvz", ",", 
             RowBox[{"Max", "[", 
              RowBox[{"Abs", "[", "\[Psi]", "]"}], "]"}], ",", 
             RowBox[{"Max", "[", 
              RowBox[{"Abs", "[", 
               RowBox[{"\[Psi]2", "-", "\[Psi]1"}], "]"}], "]"}], ",", 
             "\"\<ErrScale 1e-4\>\""}], "}"}]}], "]"}], "/.", 
         RowBox[{"b_", "\[RuleDelayed]", 
          RowBox[{
           RowBox[{"CForm", "[", "b", "]"}], "/;", 
           RowBox[{"NumericQ", "[", "b", "]"}]}]}]}]}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"$resout", "=", 
      RowBox[{"OpenAppend", "[", 
       RowBox[{"$WD", "<>", "\"\<\\\\result.dat\>\""}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"SetOptions", "[", 
      RowBox[{"$resout", ",", 
       RowBox[{"PageWidth", "\[Rule]", "\[Infinity]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"PutAppend", "[", 
      RowBox[{"putlist", ",", "$resout"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "gamma", ",", "Fr", ",", "Coriolis", ",", "shift", ",", "maxvx", ",", 
        "dvx", ",", "maxvy", ",", "dvy", ",", "maxvz", ",", "dvz", ",", 
        "mqvx", ",", "dqvx", ",", "mqvy", ",", "dqvy", ",", "mqvz", ",", 
        "dqvz", ",", "\[Psi]", ",", "\[Psi]2", ",", "\[Psi]1"}], "}"}], 
      "=."}], ";", "\[IndentingNewLine]", 
     RowBox[{"Close", "[", "$resout", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"Length", "[", "pars", "]"}]}], "}"}]}], "\[IndentingNewLine]", 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.424420539773*^9, 3.424420547672*^9}, {3.424420707809*^9,
    3.424420713389*^9}, {3.424420759052*^9, 3.4244207703389997`*^9}, 
   3.4244964698046*^9, {3.4244965404196*^9, 3.4244965907646*^9}, {
   3.424672039992*^9, 3.424672040644*^9}, {3.424685877091*^9, 
   3.424685898726*^9}, {3.4246860599700003`*^9, 3.424686089901*^9}, {
   3.4250226403256674`*^9, 3.4250226474017897`*^9}, {3.4250279418905106`*^9, 
   3.4250279748438053`*^9}, {3.425028351671485*^9, 3.4250283528966074`*^9}, {
   3.4420462538125*^9, 3.442046254796875*^9}, {3.442048608265625*^9, 
   3.44204860865625*^9}, {3.44205356553125*^9, 3.44205356628125*^9}, {
   3.4423126648125*^9, 3.442312680359375*^9}, {3.442318168484375*^9, 
   3.442318170390625*^9}, {3.4425015690810003`*^9, 3.4425015717539997`*^9}, 
   3.44250279999*^9, {3.442504533707*^9, 3.4425045570290003`*^9}, {
   3.442505197054*^9, 3.442505219318*^9}, {3.4425054209849997`*^9, 
   3.4425054225360003`*^9}, {3.442505601201*^9, 3.442505603276*^9}, {
   3.442505798258*^9, 3.442505800084*^9}, 3.4425065238059998`*^9, 
   3.4425065953199997`*^9, 3.442506774941*^9, {3.442507420257*^9, 
   3.442507436237*^9}, {3.4425080223050003`*^9, 3.442508042734*^9}, {
   3.442508112755*^9, 3.442508123516*^9}, {3.442508278344*^9, 
   3.442508279571*^9}, 3.442508681082*^9, {3.442508714898*^9, 
   3.4425087159779997`*^9}, {3.442509097354*^9, 3.442509115538*^9}, 
   3.4425091561*^9, 3.442509252835*^9, 3.442509310083*^9, 3.442509938058*^9, {
   3.442509980337*^9, 3.442509981443*^9}, {3.442510199034*^9, 
   3.442510247355*^9}, 3.4425110774*^9, {3.4425112942530003`*^9, 
   3.442511307491*^9}, {3.442511348616*^9, 3.442511349451*^9}, 
   3.4425115583310003`*^9, {3.4425116599820004`*^9, 3.442511660921*^9}, 
   3.442511717558*^9, {3.442512200615*^9, 3.4425122011619997`*^9}, {
   3.44251223701*^9, 3.442512238117*^9}, 3.442512555877*^9, 
   3.4425128875220003`*^9, 3.442514001566*^9, 3.44256471477*^9, 
   3.4426384774110003`*^9, 3.442639268377*^9, {3.442641432568*^9, 
   3.442641433236*^9}, {3.442641835356*^9, 3.4426418372720003`*^9}, 
   3.442642205058*^9, {3.442646594682*^9, 3.4426465954*^9}, {
   3.442647040994*^9, 3.4426470828789997`*^9}, {3.44264793402*^9, 
   3.442647934867*^9}, {3.4426488865690002`*^9, 3.442648886749*^9}, 
   3.442649348238*^9, {3.4426493890179996`*^9, 3.4426493899700003`*^9}, 
   3.44265027317*^9, {3.44265044346*^9, 3.4426504742539997`*^9}, {
   3.442650764574*^9, 3.4426507817720003`*^9}, {3.4428246783129997`*^9, 
   3.442824679712*^9}, {3.442825593979*^9, 3.442825611593*^9}, {
   3.442826304251*^9, 3.4428263061*^9}, 3.442826507415*^9, {
   3.4428307511610003`*^9, 3.442830752031*^9}, {3.44283107015*^9, 
   3.442831070665*^9}, {3.442831619543*^9, 3.442831620138*^9}, 
   3.4428346265030003`*^9, {3.4428351738570004`*^9, 3.4428351744309998`*^9}, {
   3.4428365579779997`*^9, 3.4428365583129997`*^9}, {3.442853403668*^9, 
   3.4428534040480003`*^9}, 3.44285369695*^9, {3.44584066375*^9, 
   3.445840664125*^9}, {3.4458407186875*^9, 3.4458407191875*^9}, {
   3.445840853984375*^9, 3.44584085521875*^9}, 3.445843044375*^9, {
   3.445843246828125*^9, 3.445843247234375*^9}, {3.445844257*^9, 
   3.445844257140625*^9}, {3.448352972538*^9, 3.4483529916005*^9}, {
   3.448694549109375*^9, 3.448694560546875*^9}, {3.448696574015625*^9, 
   3.448696575890625*^9}, {3.4486981549375*^9, 3.44869815709375*^9}, {
   3.448708845234375*^9, 3.44870884571875*^9}, {3.448709462359375*^9, 
   3.4487094626875*^9}, {3.44877703125*^9, 3.44877703184375*^9}, {
   3.4487774285625*^9, 3.44877742890625*^9}, {3.44878377809375*^9, 
   3.448783779171875*^9}, {3.448955335828125*^9, 3.448955336546875*^9}, {
   3.448955796*^9, 3.4489557978125*^9}, {3.448958334796875*^9, 
   3.448958337421875*^9}, {3.448958455515625*^9, 3.448958458015625*^9}, {
   3.44895854878125*^9, 3.448958552203125*^9}, 3.448958941453125*^9, {
   3.44895905065625*^9, 3.44895905071875*^9}, {3.4489591439375*^9, 
   3.448959146078125*^9}, {3.44895923565625*^9, 3.448959238578125*^9}, {
   3.448959301125*^9, 3.448959305109375*^9}, {3.448959353078125*^9, 
   3.448959355578125*^9}, {3.44895940759375*^9, 3.4489594100625*^9}, {
   3.448959975265625*^9, 3.448959979765625*^9}, {3.44896093871875*^9, 
   3.44896093884375*^9}, {3.44896115240625*^9, 3.448961156625*^9}, {
   3.44896353146875*^9, 3.44896353165625*^9}, {3.44896374259375*^9, 
   3.448963743015625*^9}, 3.449203645671875*^9, {3.4492118018125*^9, 
   3.4492118168125*^9}, {3.4492122759375*^9, 3.44921228315625*^9}, {
   3.44921680759375*^9, 3.449216843125*^9}, 3.449217611546875*^9, 
   3.44929218046875*^9, {3.449295434734375*^9, 3.44929546225*^9}, {
   3.449296950953125*^9, 3.4492969510625*^9}, {3.449297002046875*^9, 
   3.44929700490625*^9}, {3.449297185453125*^9, 3.449297187703125*^9}, {
   3.44929729690625*^9, 3.449297300375*^9}, {3.449297468234375*^9, 
   3.449297471046875*^9}, {3.449297553953125*^9, 3.4492975564375*^9}, {
   3.449297641921875*^9, 3.449297645671875*^9}, {3.449299009359375*^9, 
   3.449299014265625*^9}, {3.44982346515625*^9, 3.449823465421875*^9}, {
   3.449907878453125*^9, 3.449907888359375*^9}, {3.4499107680625*^9, 
   3.44991076996875*^9}, {3.4499959683125*^9, 3.4499959696875*^9}, {
   3.449996211296875*^9, 3.44999621390625*^9}, {3.449998746234375*^9, 
   3.449998750296875*^9}, {3.45007137253125*^9, 3.450071375453125*^9}, 
   3.45023583675*^9, 3.450236135609375*^9, {3.45023622740625*^9, 
   3.450236227453125*^9}, {3.45023656240625*^9, 3.4502365649375*^9}, {
   3.45041388878125*^9, 3.45041389415625*^9}, {3.45041413715625*^9, 
   3.450414137734375*^9}, {3.45041417090625*^9, 3.450414181234375*^9}, {
   3.450415025359375*^9, 3.4504150275625*^9}, {3.45041540059375*^9, 
   3.450415403328125*^9}, {3.450416030171875*^9, 3.45041603640625*^9}, {
   3.450416279515625*^9, 3.45041628184375*^9}, {3.450416388359375*^9, 
   3.450416390578125*^9}, {3.45042328396875*^9, 3.450423290859375*^9}, 
   3.458968236234375*^9, 3.4590535314375*^9, {3.459054455*^9, 
   3.45905450815625*^9}, {3.459061562640625*^9, 3.459061565078125*^9}, {
   3.4590616566875*^9, 3.459061658890625*^9}, {3.45906337284375*^9, 
   3.45906340934375*^9}, {3.4590637503125*^9, 3.459063758140625*^9}, {
   3.4590638070625*^9, 3.459063810453125*^9}, {3.45906389884375*^9, 
   3.45906390946875*^9}, {3.467008391978125*^9, 3.467008420759375*^9}, {
   3.467008476165625*^9, 3.467008476478125*^9}, {3.501441003135125*^9, 
   3.501441004916375*^9}, 3.5014827584375*^9, {3.501490815078125*^9, 
   3.50149081525*^9}, {3.50156001284375*^9, 3.501560013453125*^9}, {
   3.5015600833125*^9, 3.5015600836875*^9}, 3.501560358265625*^9, {
   3.505113046921875*^9, 3.50511306275*^9}, {3.5051180150625*^9, 
   3.50511801753125*^9}, {3.505118066703125*^9, 3.50511807196875*^9}, {
   3.505118118234375*^9, 3.505118136234375*^9}, 3.505118294375*^9, {
   3.50511879509375*^9, 3.505118796203125*^9}, {3.505118957703125*^9, 
   3.505118959125*^9}, {3.50511912059375*^9, 3.505119138953125*^9}, 
   3.505119262234375*^9, {3.50528246184375*^9, 3.50528247378125*^9}, {
   3.5052825175625*^9, 3.505282517671875*^9}, 3.505282779765625*^9, 
   3.505282938296875*^9, 3.50528306078125*^9, {3.505283458328125*^9, 
   3.50528348*^9}, {3.505291733640625*^9, 3.5052917338125*^9}, {
   3.50529179790625*^9, 3.50529180540625*^9}, 3.5052945143125*^9, 
   3.505295601703125*^9, 3.505295856453125*^9, 3.50529605653125*^9, 
   3.505296211390625*^9, {3.505296360078125*^9, 3.50529636015625*^9}, {
   3.505297056109375*^9, 3.505297056171875*^9}, 3.505297159703125*^9, 
   3.50529720196875*^9, 3.50529725765625*^9, {3.5052994333125*^9, 
   3.50529943346875*^9}, {3.505299467484375*^9, 3.505299512*^9}, 
   3.505299577765625*^9, 3.505299613703125*^9, 3.50529966634375*^9, 
   3.505299730640625*^9, {3.50535511421875*^9, 3.50535511865625*^9}, 
   3.50535541440625*^9, {3.505355467078125*^9, 3.505355477578125*^9}, {
   3.505357024921875*^9, 3.505357044921875*^9}, {3.5053570835*^9, 
   3.505357219109375*^9}, {3.505357338765625*^9, 3.5053575203125*^9}, 
   3.505357602359375*^9, {3.50535810834375*^9, 3.50535813946875*^9}, {
   3.505358394796875*^9, 3.5053584214375*^9}, {3.505358804375*^9, 
   3.50535881421875*^9}, {3.505358952203125*^9, 3.505358973203125*^9}, {
   3.505359078140625*^9, 3.505359122578125*^9}, {3.505359164984375*^9, 
   3.505359165109375*^9}, {3.505359248328125*^9, 3.505359249859375*^9}, {
   3.505359452953125*^9, 3.505359458109375*^9}, {3.505359490375*^9, 
   3.505359526125*^9}, {3.505400702890625*^9, 3.505400704640625*^9}, {
   3.505406893046875*^9, 3.5054068966875*^9}, {3.505442598890625*^9, 
   3.50544260890625*^9}, {3.5054451786875*^9, 3.505445186890625*^9}, {
   3.50544936884375*^9, 3.50544938628125*^9}, {3.5054497074375*^9, 
   3.505449731984375*^9}, {3.505449792859375*^9, 3.50544984690625*^9}, {
   3.505449993375*^9, 3.50545000253125*^9}, {3.505450057203125*^9, 
   3.50545006328125*^9}, {3.505450154*^9, 3.505450159546875*^9}, {
   3.50545035234375*^9, 3.505450357078125*^9}, {3.505450426*^9, 
   3.505450426859375*^9}, {3.50545066671875*^9, 3.505450670546875*^9}, {
   3.505450710234375*^9, 3.505450711984375*^9}, {3.505450840015625*^9, 
   3.505450840671875*^9}, {3.50545765009375*^9, 3.50545769615625*^9}, {
   3.505457950609375*^9, 3.505457950859375*^9}, {3.505458372125*^9, 
   3.5054583735*^9}, {3.505458766984375*^9, 3.505458771140625*^9}, {
   3.5054605455*^9, 3.5054605469375*^9}, {3.50546065390625*^9, 
   3.505460657*^9}, {3.50546209325*^9, 3.50546209721875*^9}, {
   3.505462560859375*^9, 3.5054625669375*^9}, {3.5054634454375*^9, 
   3.505463473203125*^9}, 3.505463539421875*^9, {3.505465609046875*^9, 
   3.505465621953125*^9}, {3.50546604915625*^9, 3.505466049203125*^9}, {
   3.505466238796875*^9, 3.50546624878125*^9}, {3.505468282171875*^9, 
   3.505468284140625*^9}, {3.505468323734375*^9, 3.505468332953125*^9}, {
   3.50546837353125*^9, 3.505468373609375*^9}, {3.50614766903275*^9, 
   3.50614767034525*^9}, 3.50614771290775*^9, {3.506148037392125*^9, 
   3.5061480480015*^9}, {3.5061484353765*^9, 3.506148446954625*^9}, {
   3.506149208798375*^9, 3.50614921097025*^9}, {3.506150064985875*^9, 
   3.506150069735875*^9}, {3.50632120590775*^9, 3.506321207439*^9}, {
   3.506664842142125*^9, 3.5066648721265*^9}, {3.5066649803765*^9, 
   3.50666498140775*^9}, {3.506665037314*^9, 3.50666513084525*^9}, {
   3.5069207033805*^9, 3.506920717271125*^9}, {3.506920807083625*^9, 
   3.50692081241175*^9}, {3.50692457416175*^9, 3.506924576396125*^9}, {
   3.50692468828675*^9, 3.506924699708625*^9}, {3.506924762708625*^9, 
   3.506924845271125*^9}, {3.50692504816175*^9, 3.5069250501305*^9}, {
   3.50726579209925*^9, 3.50726579341175*^9}, {3.50726597703675*^9, 
   3.507265978052375*^9}, {3.507266184583625*^9, 3.507266209677375*^9}, 
   3.507266330568*^9, 3.50726636109925*^9, {3.507266601958625*^9, 
   3.507266603239875*^9}, {3.50726668303675*^9, 3.507266684396125*^9}, {
   3.50726675103675*^9, 3.50726676891175*^9}, {3.507267430958625*^9, 
   3.507267431083625*^9}, {3.5072681642555*^9, 3.5072681852555*^9}, {
   3.5072687438805*^9, 3.50726875234925*^9}, {3.507269760271125*^9, 
   3.507269766396125*^9}, 3.507269807177375*^9, {3.5072706378805*^9, 
   3.507270644568*^9}, {3.50728151903675*^9, 3.507281525083625*^9}, 
   3.5072816558805*^9, {3.50728224541175*^9, 3.507282252364875*^9}, 
   3.50728237853675*^9, {3.50728290734925*^9, 3.507282908568*^9}, 
   3.536049312431506*^9, {3.536049357744006*^9, 3.536049433244006*^9}, {
   3.5360494901658816`*^9, 3.536049578556506*^9}, 3.536049772212756*^9, {
   3.5360499023846316`*^9, 3.5360499034158816`*^9}, 3.5360499380408816`*^9, {
   3.536050075337756*^9, 3.5360500911971316`*^9}, {3.538384632640625*^9, 
   3.5383846836875*^9}, {3.53838471803125*^9, 3.53838473909375*^9}, 
   3.538384857234375*^9, {3.538385025984375*^9, 3.53838503109375*^9}, {
   3.5383851689375*^9, 3.5383851760625*^9}, {3.538385271234375*^9, 
   3.538385333921875*^9}, 3.53838547653125*^9, 3.53838555525*^9, {
   3.538385846125*^9, 3.538385859828125*^9}, {3.538385894515625*^9, 
   3.538385896859375*^9}, {3.538385930515625*^9, 3.538385961*^9}, 
   3.5383859910625*^9, {3.53838602515625*^9, 3.538386045453125*^9}, {
   3.53838616103125*^9, 3.538386161296875*^9}, {3.53838621303125*^9, 
   3.5383862245*^9}, {3.538386368375*^9, 3.5383863688125*^9}, {
   3.53838657903125*^9, 3.538386595921875*^9}, {3.5387380713125*^9, 
   3.5387380926875*^9}, {3.53873829765625*^9, 3.53873829809375*^9}, {
   3.539334962817375*^9, 3.53933496405175*^9}, {3.5998130324873123`*^9, 
   3.599813065052742*^9}, {3.599813729537117*^9, 3.5998137296904373`*^9}, {
   3.5998160101503983`*^9, 3.599816018771492*^9}, {3.5998949365400467`*^9, 
   3.599894955188484*^9}, 3.6036179341128206`*^9, {3.604213491459961*^9, 
   3.6042135146728516`*^9}, {3.6042139439541016`*^9, 
   3.6042140410117188`*^9}, {3.6042140955683594`*^9, 3.604214127256836*^9}, {
   3.6042148426845703`*^9, 3.6042148685810547`*^9}, {3.604215026104492*^9, 
   3.6042150485507812`*^9}, {3.604215332701172*^9, 3.6042153549189453`*^9}, {
   3.604299514930664*^9, 3.6042995150351562`*^9}, {3.604299603953125*^9, 
   3.604299609192383*^9}, {3.6042996486992188`*^9, 3.6042996517158203`*^9}, 
   3.605260594581328*^9, {3.605260641756133*^9, 3.6052606474377737`*^9}, {
   3.6052607039387503`*^9, 3.605260704365508*^9}, {3.6052607351369925`*^9, 
   3.6052607893752737`*^9}, {3.605260855466094*^9, 3.605260902464141*^9}, {
   3.605261026493438*^9, 3.605261064921172*^9}, {3.605261116849883*^9, 
   3.6052611388869925`*^9}, {3.605261428195586*^9, 3.605261433612578*^9}, 
   3.6052623950383596`*^9, {3.6052632217395315`*^9, 3.605263227149688*^9}, {
   3.605263685577422*^9, 3.6052636856887503`*^9}, {3.623483844267578*^9, 
   3.623483847326172*^9}, 3.6234839614814453`*^9, {3.6234840821445312`*^9, 
   3.623484086267578*^9}, {3.6234847819472656`*^9, 3.623484785084961*^9}, {
   3.6234848228085938`*^9, 3.6234848444726562`*^9}, {3.6234870516328125`*^9, 
   3.6234870679833984`*^9}, {3.6235552766992188`*^9, 
   3.6235552771992188`*^9}, {3.623652368368164*^9, 3.623652377520508*^9}, {
   3.6236543717929688`*^9, 3.623654372586914*^9}, {3.624259724213867*^9, 
   3.6242597358964844`*^9}, {3.6242598083779297`*^9, 3.624259817979492*^9}, {
   3.6246883643496094`*^9, 3.6246883791845703`*^9}, {3.624697689006836*^9, 
   3.6246976890722656`*^9}, {3.6246977864726562`*^9, 
   3.6246977955214844`*^9}, {3.6248644488808594`*^9, 3.624864463741211*^9}, {
   3.624866482055664*^9, 3.6248664958242188`*^9}, {3.625906809404297*^9, 
   3.625906841607422*^9}, 3.626072846783203*^9, {3.6260742931132812`*^9, 
   3.6260742966572266`*^9}, {3.6260749544257812`*^9, 
   3.6260749642773438`*^9}, {3.6260769987158203`*^9, 3.626076999602539*^9}, {
   3.626079964850586*^9, 3.6260799871484375`*^9}, {3.6266125129785156`*^9, 
   3.6266125290722656`*^9}, {3.6266224870039062`*^9, 3.626622497448242*^9}, 
   3.626622749901367*^9, {3.626625101629883*^9, 3.6266251036503906`*^9}, {
   3.6266258475966797`*^9, 3.626625849236328*^9}, {3.626626442439453*^9, 
   3.626626451526367*^9}, {3.626626775711914*^9, 3.6266267773945312`*^9}, {
   3.6266271591953125`*^9, 3.626627161263672*^9}, {3.6266280307841797`*^9, 
   3.6266280351914062`*^9}, 3.626629667794922*^9, {3.626630132963867*^9, 
   3.6266301436376953`*^9}, {3.626631402508789*^9, 3.626631404388672*^9}, {
   3.6266326508710938`*^9, 3.6266326869052734`*^9}, {3.6266328176591797`*^9, 
   3.6266328181152344`*^9}, 3.6266332096240234`*^9, 3.6266658182753906`*^9, {
   3.6266680730029297`*^9, 3.6266680743779297`*^9}, {3.626670824567383*^9, 
   3.626670831864258*^9}}]
}, Closed]]
},
WindowSize->{1006, 662},
WindowMargins->{{Automatic, 15}, {Automatic, 16}},
ShowSelection->True,
Magnification->1.1000001430511475`,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 120, 2, 32, "Input"],
Cell[668, 24, 144, 3, 32, "Input"],
Cell[815, 29, 1284, 18, 32, "Input"],
Cell[2102, 49, 199, 3, 32, "Input"],
Cell[2304, 54, 1203, 25, 73, "Input"],
Cell[CellGroupData[{
Cell[3532, 83, 205, 2, 99, "Title"],
Cell[3740, 87, 29641, 640, 1320, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[33418, 732, 205, 2, 70, "Title"],
Cell[33626, 736, 29870, 624, 1252, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[63533, 1365, 194, 2, 70, "Title"],
Cell[63730, 1369, 34058, 685, 1352, "Input"]
}, Closed]]
}
]
*)

(* End of internal cache information *)
